Main network model diagnostics - balanced statistics

This file shows diagnostics for main network models fit using balanced racial/ethnic mixing matrices and degree terms adjusted to correspond to the balanced mixing matrices. In this file, we fit a series of nested models by adding one term at a time to examine changes to model estimates, MCMC diagnostics, and network diagnostics.

Load packages and model fits

rm(list = ls())
suppressMessages(library("EpiModelHIV"))
library("latticeExtra")
## Loading required package: lattice
## Loading required package: RColorBrewer
library("knitr")
library("kableExtra")

load(file = "/homes/dpwhite/R/GitHub Repos/WHAMP/Model fits and simulations/Fit tests and debugging/est/fit.m.buildup.bal.rda")

Model terms and control settings

Model terms and target statistics
Terms Model 1 Model 2 Model 3 Model 4 Model 5 Model 6 Model 7 Model 8
edges 2240.5 2240.5 2240.5 2240.5 2240.5 2240.5 2240.5 2240.5
nodefactor.deg.pers.1 NA NA NA 493.0 493.0 493.0 493.0 493.0
nodefactor.deg.pers.2 NA NA NA 603.0 603.0 603.0 603.0 603.0
nodefactor.race..wa.B NA 213.8 213.8 213.8 213.8 213.8 213.8 213.8
nodefactor.race..wa.H NA 587.8 587.8 587.8 587.8 587.8 587.8 587.8
nodefactor.region.EW NA NA NA NA 445.6 445.6 445.6 445.6
nodefactor.region.OW NA NA NA NA 1278.1 1278.1 1278.1 1278.1
nodematch.race..wa.B NA NA 31.2 31.2 31.2 31.2 31.2 31.2
nodematch.race..wa.H NA NA 123.3 123.3 123.3 123.3 123.3 123.3
nodematch.race..wa.O NA NA 1638.9 1638.9 1638.9 1638.9 1638.9 1638.9
absdiff.sqrt.age NA NA NA NA NA 1206.3 1206.3 1206.3
degrange 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
nodematch.role.class.I -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf
nodematch.role.class.R -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf
mix.region.EW.KC NA NA NA NA NA NA -Inf NA
mix.region.EW.OW NA NA NA NA NA NA -Inf NA
mix.region.KC.OW NA NA NA NA NA NA -Inf NA
nodematch.region NA NA NA NA NA NA NA 2016.5

The control settings for these models are:

set.control.ergm = control.ergm(MCMC.interval = 1e+5,
                                MCMC.samplesize = 7500,
                                MCMC.burnin = 1e+6,
                                MPLE.max.dyad.types = 1e+7,
                                init.method = "zeros",
                                MCMLE.maxit = 400,
                                parallel = np/2, 
                                parallel.type="PSOCK"))

MCMC diagnostics

Model 1

## Sample statistics summary:
## 
## Iterations = 1e+06:750900000
## Thinning interval = 1e+05 
## Number of chains = 4 
## Sample size per chain = 7500 
## 
## 1. Empirical mean and standard deviation for each variable,
##    plus standard error of the mean:
## 
##           Mean             SD       Naive SE Time-series SE 
##        0.04827       29.26193        0.16894        0.17268 
## 
## 2. Quantiles for each variable:
## 
##  2.5%   25%   50%   75% 97.5% 
## -57.5 -19.5   0.5  19.5  57.5 
## 
## 
## Sample statistics cross-correlations:
##       edges
## edges     1
## 
## Sample statistics auto-correlation:
## Chain 1 
##                  edges
## Lag 0      1.000000000
## Lag 1e+05  0.000812571
## Lag 2e+05 -0.014008792
## Lag 3e+05 -0.013181798
## Lag 4e+05 -0.024329088
## Lag 5e+05 -0.004715004
## Chain 2 
##                   edges
## Lag 0      1.0000000000
## Lag 1e+05  0.0113871359
## Lag 2e+05  0.0051766637
## Lag 3e+05  0.0095815335
## Lag 4e+05  0.0018761518
## Lag 5e+05 -0.0007773829
## Chain 3 
##                  edges
## Lag 0      1.000000000
## Lag 1e+05  0.007762176
## Lag 2e+05  0.028607313
## Lag 3e+05 -0.007422592
## Lag 4e+05 -0.012009726
## Lag 5e+05 -0.016179124
## Chain 4 
##                edges
## Lag 0     1.00000000
## Lag 1e+05 0.01722746
## Lag 2e+05 0.02051297
## Lag 3e+05 0.00415593
## Lag 4e+05 0.01637713
## Lag 5e+05 0.01789402
## 
## Sample statistics burn-in diagnostic (Geweke):
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##  edges 
## 0.8729 
## 
## Individual P-values (lower = worse):
##     edges 
## 0.3827018 
## Joint P-value (lower = worse):  0.3787341 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##  edges 
## -1.239 
## 
## Individual P-values (lower = worse):
##     edges 
## 0.2151753 
## Joint P-value (lower = worse):  0.2518375 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##  edges 
## 0.7333 
## 
## Individual P-values (lower = worse):
##     edges 
## 0.4633587 
## Joint P-value (lower = worse):  0.4428171 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##   edges 
## -0.4826 
## 
## Individual P-values (lower = worse):
##     edges 
## 0.6293818 
## Joint P-value (lower = worse):  0.6262199 .
## Warning in formals(fun): argument is not a function

## 
## MCMC diagnostics shown here are from the last round of simulation, prior to computation of final parameter estimates. Because the final estimates are refinements of those used for this simulation run, these diagnostics may understate model performance. To directly assess the performance of the final model on in-model statistics, please use the GOF command: gof(ergmFitObject, GOF=~model).

Model 2

## Sample statistics summary:
## 
## Iterations = 1e+06:750900000
## Thinning interval = 1e+05 
## Number of chains = 4 
## Sample size per chain = 7500 
## 
## 1. Empirical mean and standard deviation for each variable,
##    plus standard error of the mean:
## 
##                           Mean    SD Naive SE Time-series SE
## edges                  0.16140 28.89  0.16677        0.16761
## nodefactor.race..wa.B -0.35277 11.93  0.06889        0.06860
## nodefactor.race..wa.H -0.03357 16.51  0.09535        0.09635
## 
## 2. Quantiles for each variable:
## 
##                         2.5%     25%    50%    75% 97.5%
## edges                 -57.50 -19.500 0.5000 19.500 56.50
## nodefactor.race..wa.B -23.83  -8.834 0.1664  7.166 23.17
## nodefactor.race..wa.H -31.84 -10.844 0.1560 11.156 32.16
## 
## 
## Sample statistics cross-correlations:
##                           edges nodefactor.race..wa.B
## edges                 1.0000000            0.27886691
## nodefactor.race..wa.B 0.2788669            1.00000000
## nodefactor.race..wa.H 0.3788393            0.02229181
##                       nodefactor.race..wa.H
## edges                            0.37883930
## nodefactor.race..wa.B            0.02229181
## nodefactor.race..wa.H            1.00000000
## 
## Sample statistics auto-correlation:
## Chain 1 
##                   edges nodefactor.race..wa.B nodefactor.race..wa.H
## Lag 0      1.0000000000           1.000000000           1.000000000
## Lag 1e+05  0.0195321964           0.006087381           0.013009692
## Lag 2e+05 -0.0053867865           0.020550132          -0.004257186
## Lag 3e+05  0.0007608276           0.001542963          -0.005080063
## Lag 4e+05  0.0139952676          -0.008069587          -0.002986683
## Lag 5e+05 -0.0049831999          -0.011553817          -0.011021465
## Chain 2 
##                  edges nodefactor.race..wa.B nodefactor.race..wa.H
## Lag 0      1.000000000          1.0000000000           1.000000000
## Lag 1e+05  0.006972719          0.0001342158           0.022476011
## Lag 2e+05 -0.002133960          0.0056226991          -0.007535108
## Lag 3e+05 -0.005716569         -0.0223128768           0.009045639
## Lag 4e+05  0.014376243         -0.0213538759           0.007926069
## Lag 5e+05 -0.002807807          0.0074446690           0.001725296
## Chain 3 
##                  edges nodefactor.race..wa.B nodefactor.race..wa.H
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05 -0.006531002          -0.016808133          -0.003807935
## Lag 2e+05  0.008894848           0.001996580          -0.012241177
## Lag 3e+05 -0.005820613           0.007095523          -0.001577152
## Lag 4e+05 -0.008469688          -0.003821740          -0.022402398
## Lag 5e+05 -0.007095829           0.015782163          -0.005541702
## Chain 4 
##                   edges nodefactor.race..wa.B nodefactor.race..wa.H
## Lag 0      1.000000e+00           1.000000000           1.000000000
## Lag 1e+05  2.515585e-03           0.008988164           0.019596458
## Lag 2e+05  1.723743e-02          -0.012986270          -0.015809080
## Lag 3e+05 -1.198866e-02           0.010389819           0.006600693
## Lag 4e+05 -1.834574e-05           0.015947977          -0.010012109
## Lag 5e+05 -7.682528e-03           0.003210243          -0.004021384
## 
## Sample statistics burn-in diagnostic (Geweke):
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##                1.2078               -1.3891                0.1263 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##             0.2271064             0.1647872             0.8995327 
## Joint P-value (lower = worse):  0.1626046 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##               -1.0034               -1.8141               -0.9348 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##            0.31565589            0.06965786            0.34990654 
## Joint P-value (lower = worse):  0.2638284 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##               -0.5860               -1.6830               -0.7085 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##            0.55789177            0.09238176            0.47865467 
## Joint P-value (lower = worse):  0.3624487 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##               -0.3623                1.1003                0.2287 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##             0.7171220             0.2711948             0.8191352 
## Joint P-value (lower = worse):  0.5859225 .
## Warning in formals(fun): argument is not a function

## 
## MCMC diagnostics shown here are from the last round of simulation, prior to computation of final parameter estimates. Because the final estimates are refinements of those used for this simulation run, these diagnostics may understate model performance. To directly assess the performance of the final model on in-model statistics, please use the GOF command: gof(ergmFitObject, GOF=~model).

Model 3

## Sample statistics summary:
## 
## Iterations = 1e+06:750900000
## Thinning interval = 1e+05 
## Number of chains = 4 
## Sample size per chain = 7500 
## 
## 1. Empirical mean and standard deviation for each variable,
##    plus standard error of the mean:
## 
##                           Mean     SD Naive SE Time-series SE
## edges                 -0.04597 29.034  0.16763        0.16763
## nodefactor.race..wa.B  0.10420 12.699  0.07332        0.07490
## nodefactor.race..wa.H  0.29213 17.610  0.10167        0.10904
## nodematch.race..wa.B   0.11116  4.981  0.02876        0.03089
## nodematch.race..wa.H   0.27055  8.757  0.05056        0.06556
## nodematch.race..wa.O  -0.22163 26.416  0.15251        0.15232
## 
## 2. Quantiles for each variable:
## 
##                          2.5%     25%      50%    75%  97.5%
## edges                 -57.500 -19.500 -0.50000 19.500 56.500
## nodefactor.race..wa.B -24.834  -8.834  0.16640  8.166 25.166
## nodefactor.race..wa.H -33.844 -11.844  0.15600 12.156 35.156
## nodematch.race..wa.B   -9.177  -3.177 -0.17694  3.823  9.823
## nodematch.race..wa.H  -16.300  -5.300  0.69972  5.700 17.700
## nodematch.race..wa.O  -51.946 -17.946  0.05397 17.054 52.054
## 
## 
## Sample statistics cross-correlations:
##                           edges nodefactor.race..wa.B
## edges                 1.0000000            0.26341218
## nodefactor.race..wa.B 0.2634122            1.00000000
## nodefactor.race..wa.H 0.3600299            0.04743928
## nodematch.race..wa.B  0.1125557            0.57542603
## nodematch.race..wa.H  0.1304341           -0.04962059
## nodematch.race..wa.O  0.8191997           -0.05277563
##                       nodefactor.race..wa.H nodematch.race..wa.B
## edges                            0.36002992           0.11255571
## nodefactor.race..wa.B            0.04743928           0.57542603
## nodefactor.race..wa.H            1.00000000          -0.01750918
## nodematch.race..wa.B            -0.01750918           1.00000000
## nodematch.race..wa.H             0.59476949           0.01125704
## nodematch.race..wa.O            -0.06283077           0.02212869
##                       nodematch.race..wa.H nodematch.race..wa.O
## edges                           0.13043412           0.81919975
## nodefactor.race..wa.B          -0.04962059          -0.05277563
## nodefactor.race..wa.H           0.59476949          -0.06283077
## nodematch.race..wa.B            0.01125704           0.02212869
## nodematch.race..wa.H            1.00000000           0.07138299
## nodematch.race..wa.O            0.07138299           1.00000000
## 
## Sample statistics auto-correlation:
## Chain 1 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H
## Lag 0     1.000000000           1.000000000           1.000000000
## Lag 1e+05 0.011208103           0.029789555           0.054696518
## Lag 2e+05 0.005052589          -0.010114664           0.014285503
## Lag 3e+05 0.002312177          -0.007342714          -0.023233439
## Lag 4e+05 0.004227302           0.014950665          -0.020507712
## Lag 5e+05 0.003788233           0.003250517          -0.002222634
##           nodematch.race..wa.B nodematch.race..wa.H nodematch.race..wa.O
## Lag 0              1.000000000          1.000000000         1.0000000000
## Lag 1e+05          0.052004416          0.242352007         0.0116599475
## Lag 2e+05         -0.011522050          0.064701631        -0.0033268290
## Lag 3e+05         -0.003198173          0.004514173         0.0175037043
## Lag 4e+05          0.001351352         -0.006320866         0.0021363092
## Lag 5e+05         -0.021783436         -0.013148267         0.0000869757
## Chain 2 
##                  edges nodefactor.race..wa.B nodefactor.race..wa.H
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05  0.008315285           0.023384221           0.076628361
## Lag 2e+05  0.005919698          -0.009589067           0.030624057
## Lag 3e+05 -0.004743854          -0.007090496          -0.005666997
## Lag 4e+05 -0.020476812          -0.003349672           0.012568118
## Lag 5e+05  0.009085159           0.013542092          -0.002743474
##           nodematch.race..wa.B nodematch.race..wa.H nodematch.race..wa.O
## Lag 0              1.000000000         1.0000000000          1.000000000
## Lag 1e+05          0.053192246         0.2640043133         -0.007635212
## Lag 2e+05          0.010216815         0.0819883136          0.011403395
## Lag 3e+05          0.018601429         0.0120562369          0.005709449
## Lag 4e+05          0.005380062         0.0234939229         -0.011390689
## Lag 5e+05          0.027208215         0.0005846989          0.013501395
## Chain 3 
##                   edges nodefactor.race..wa.B nodefactor.race..wa.H
## Lag 0      1.0000000000           1.000000000           1.000000000
## Lag 1e+05 -0.0025697584           0.032638850           0.074429674
## Lag 2e+05 -0.0003722483           0.011145369           0.010129884
## Lag 3e+05  0.0005192685           0.002223597           0.021366519
## Lag 4e+05  0.0112559960          -0.007507226          -0.006198121
## Lag 5e+05  0.0014807265           0.004903909          -0.001255750
##           nodematch.race..wa.B nodematch.race..wa.H nodematch.race..wa.O
## Lag 0              1.000000000          1.000000000          1.000000000
## Lag 1e+05          0.076422422          0.247745845         -0.016008765
## Lag 2e+05          0.013355573          0.075227808          0.003136806
## Lag 3e+05          0.022249190          0.018094273          0.004314354
## Lag 4e+05         -0.006144743         -0.002030739          0.012968243
## Lag 5e+05         -0.007460851          0.001254804         -0.007126456
## Chain 4 
##                  edges nodefactor.race..wa.B nodefactor.race..wa.H
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05  0.005985276           0.007950551           0.075829176
## Lag 2e+05  0.015745851           0.005858133           0.012460052
## Lag 3e+05 -0.003117628           0.006828756           0.002508038
## Lag 4e+05 -0.009449592          -0.011924740           0.004849127
## Lag 5e+05  0.002574178          -0.014319588          -0.005092043
##           nodematch.race..wa.B nodematch.race..wa.H nodematch.race..wa.O
## Lag 0              1.000000000          1.000000000          1.000000000
## Lag 1e+05          0.047420692          0.255995352         -0.007210961
## Lag 2e+05          0.006274315          0.046861111          0.020456257
## Lag 3e+05          0.015714298          0.004293817         -0.018927599
## Lag 4e+05         -0.011147433          0.016919669          0.001546832
## Lag 5e+05         -0.008059116          0.010642682          0.011775199
## 
## Sample statistics burn-in diagnostic (Geweke):
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##                1.9958                1.8304                0.4029 
##  nodematch.race..wa.B  nodematch.race..wa.H  nodematch.race..wa.O 
##                0.1850               -0.6769                1.2596 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##            0.04595596            0.06719454            0.68699311 
##  nodematch.race..wa.B  nodematch.race..wa.H  nodematch.race..wa.O 
##            0.85320383            0.49849653            0.20782123 
## Joint P-value (lower = worse):  0.2286995 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##               -0.4395                0.7850               -1.8789 
##  nodematch.race..wa.B  nodematch.race..wa.H  nodematch.race..wa.O 
##                0.9023               -2.3143                0.3429 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##            0.66030447            0.43247516            0.06026413 
##  nodematch.race..wa.B  nodematch.race..wa.H  nodematch.race..wa.O 
##            0.36691689            0.02064973            0.73164045 
## Joint P-value (lower = worse):  0.1785632 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##                0.2101               -1.5322                0.7221 
##  nodematch.race..wa.B  nodematch.race..wa.H  nodematch.race..wa.O 
##               -0.2804                0.1064                0.2287 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##             0.8336191             0.1254656             0.4702420 
##  nodematch.race..wa.B  nodematch.race..wa.H  nodematch.race..wa.O 
##             0.7791342             0.9152534             0.8191114 
## Joint P-value (lower = worse):  0.6427681 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##               0.84547              -1.86946               0.08127 
##  nodematch.race..wa.B  nodematch.race..wa.H  nodematch.race..wa.O 
##              -1.43537              -0.15651               1.75705 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.race..wa.B nodefactor.race..wa.H 
##            0.39785073            0.06155932            0.93522395 
##  nodematch.race..wa.B  nodematch.race..wa.H  nodematch.race..wa.O 
##            0.15118048            0.87563043            0.07890863 
## Joint P-value (lower = worse):  0.1701983 .
## Warning in formals(fun): argument is not a function

## 
## MCMC diagnostics shown here are from the last round of simulation, prior to computation of final parameter estimates. Because the final estimates are refinements of those used for this simulation run, these diagnostics may understate model performance. To directly assess the performance of the final model on in-model statistics, please use the GOF command: gof(ergmFitObject, GOF=~model).

Model 4

## Sample statistics summary:
## 
## Iterations = 1e+06:750900000
## Thinning interval = 1e+05 
## Number of chains = 4 
## Sample size per chain = 7500 
## 
## 1. Empirical mean and standard deviation for each variable,
##    plus standard error of the mean:
## 
##                          Mean     SD Naive SE Time-series SE
## edges                  1.1409 29.135  0.16821        0.16821
## nodefactor.deg.pers.1  0.3364 17.719  0.10230        0.10019
## nodefactor.deg.pers.2 -0.5685 18.867  0.10893        0.11298
## nodefactor.race..wa.B  0.3244 12.728  0.07348        0.07535
## nodefactor.race..wa.H  1.4040 17.540  0.10127        0.10912
## nodematch.race..wa.B   0.1903  5.069  0.02926        0.03154
## nodematch.race..wa.H   0.5690  8.759  0.05057        0.06619
## nodematch.race..wa.O  -0.1223 26.512  0.15307        0.15307
## 
## 2. Quantiles for each variable:
## 
##                          2.5%     25%        50%    75% 97.5%
## edges                 -55.500 -18.500  1.500e+00 20.500 58.50
## nodefactor.deg.pers.1 -34.000 -11.000  5.684e-14 12.000 35.00
## nodefactor.deg.pers.2 -38.000 -13.000 -1.000e+00 12.000 36.00
## nodefactor.race..wa.B -24.834  -7.834  1.664e-01  9.166 25.17
## nodefactor.race..wa.H -32.844  -9.844  1.156e+00 13.156 35.16
## nodematch.race..wa.B   -9.177  -3.177 -1.769e-01  3.823 10.82
## nodematch.race..wa.H  -16.300  -5.300  6.997e-01  6.700 17.70
## nodematch.race..wa.O  -50.946 -17.946  5.397e-02 17.054 52.05
## 
## 
## Sample statistics cross-correlations:
##                           edges nodefactor.deg.pers.1
## edges                 1.0000000            0.40132066
## nodefactor.deg.pers.1 0.4013207            1.00000000
## nodefactor.deg.pers.2 0.4227455            0.05234482
## nodefactor.race..wa.B 0.2702926            0.10979860
## nodefactor.race..wa.H 0.3606214            0.15682145
## nodematch.race..wa.B  0.1157809            0.04630690
## nodematch.race..wa.H  0.1372652            0.05534839
## nodematch.race..wa.O  0.8226218            0.32384511
##                       nodefactor.deg.pers.2 nodefactor.race..wa.B
## edges                            0.42274555            0.27029260
## nodefactor.deg.pers.1            0.05234482            0.10979860
## nodefactor.deg.pers.2            1.00000000            0.13133040
## nodefactor.race..wa.B            0.13133040            1.00000000
## nodefactor.race..wa.H            0.17403882            0.04682189
## nodematch.race..wa.B             0.06432286            0.57418114
## nodematch.race..wa.H             0.07169058           -0.03696470
## nodematch.race..wa.O             0.33529678           -0.03963260
##                       nodefactor.race..wa.H nodematch.race..wa.B
## edges                            0.36062137           0.11578086
## nodefactor.deg.pers.1            0.15682145           0.04630690
## nodefactor.deg.pers.2            0.17403882           0.06432286
## nodefactor.race..wa.B            0.04682189           0.57418114
## nodefactor.race..wa.H            1.00000000          -0.01569720
## nodematch.race..wa.B            -0.01569720           1.00000000
## nodematch.race..wa.H             0.59798072           0.02798280
## nodematch.race..wa.O            -0.05944949           0.02943175
##                       nodematch.race..wa.H nodematch.race..wa.O
## edges                           0.13726524           0.82262182
## nodefactor.deg.pers.1           0.05534839           0.32384511
## nodefactor.deg.pers.2           0.07169058           0.33529678
## nodefactor.race..wa.B          -0.03696470          -0.03963260
## nodefactor.race..wa.H           0.59798072          -0.05944949
## nodematch.race..wa.B            0.02798280           0.02943175
## nodematch.race..wa.H            1.00000000           0.07249846
## nodematch.race..wa.O            0.07249846           1.00000000
## 
## Sample statistics auto-correlation:
## Chain 1 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000          1.0000000000           1.000000000
## Lag 1e+05  0.002745247         -0.0005698549           0.005736073
## Lag 2e+05 -0.016436802         -0.0038128123          -0.018760251
## Lag 3e+05  0.011433790          0.0013849737          -0.014862439
## Lag 4e+05  0.008482136         -0.0195622978          -0.011919800
## Lag 5e+05  0.019341633          0.0016227486           0.001014457
##           nodefactor.race..wa.B nodefactor.race..wa.H nodematch.race..wa.B
## Lag 0               1.000000000           1.000000000          1.000000000
## Lag 1e+05           0.020136029           0.071397656          0.067273758
## Lag 2e+05          -0.009847821           0.033410295          0.003976353
## Lag 3e+05          -0.000639021           0.003845062         -0.003821269
## Lag 4e+05           0.009923656          -0.007287995         -0.001324894
## Lag 5e+05          -0.011802051           0.001544721          0.003124793
##           nodematch.race..wa.H nodematch.race..wa.O
## Lag 0               1.00000000          1.000000000
## Lag 1e+05           0.26893548          0.003606402
## Lag 2e+05           0.07623866         -0.007742246
## Lag 3e+05           0.02484582          0.008752092
## Lag 4e+05           0.02287966         -0.007704242
## Lag 5e+05           0.02220601          0.017808442
## Chain 2 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000          1.0000000000
## Lag 1e+05  0.003163241           0.007325812          0.0149381762
## Lag 2e+05 -0.001253138          -0.002291310         -0.0016429226
## Lag 3e+05  0.005139615           0.008733096          0.0021712699
## Lag 4e+05 -0.020638224          -0.024516227         -0.0003938711
## Lag 5e+05  0.008635419          -0.005315609          0.0368029299
##           nodefactor.race..wa.B nodefactor.race..wa.H nodematch.race..wa.B
## Lag 0              1.0000000000          1.0000000000          1.000000000
## Lag 1e+05          0.0317485820          0.0725305787          0.081356012
## Lag 2e+05         -0.0066859368          0.0057115100          0.017594969
## Lag 3e+05         -0.0189778087         -0.0001302909         -0.007434138
## Lag 4e+05          0.0002408441         -0.0033021870          0.002455499
## Lag 5e+05         -0.0230510572          0.0015200933         -0.004694521
##           nodematch.race..wa.H nodematch.race..wa.O
## Lag 0              1.000000000         1.0000000000
## Lag 1e+05          0.253929656         0.0007777483
## Lag 2e+05          0.050732118        -0.0048771621
## Lag 3e+05          0.004803744         0.0094859891
## Lag 4e+05          0.004749943        -0.0110967666
## Lag 5e+05         -0.005699775         0.0083651479
## Chain 3 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05  0.003360710           0.008422362           0.006027530
## Lag 2e+05  0.006978230           0.005955221          -0.015985002
## Lag 3e+05 -0.024733098           0.001614016          -0.014265209
## Lag 4e+05  0.005719774          -0.005355170           0.025552569
## Lag 5e+05  0.016981669           0.005586289           0.007752374
##           nodefactor.race..wa.B nodefactor.race..wa.H nodematch.race..wa.B
## Lag 0               1.000000000           1.000000000         1.0000000000
## Lag 1e+05           0.019272486           0.043146518         0.0772420620
## Lag 2e+05           0.003247312           0.017050074         0.0206703658
## Lag 3e+05          -0.024403018           0.010747317        -0.0020094039
## Lag 4e+05           0.001187666          -0.004010833         0.0004618976
## Lag 5e+05          -0.011393580           0.012837439        -0.0035139164
##           nodematch.race..wa.H nodematch.race..wa.O
## Lag 0              1.000000000          1.000000000
## Lag 1e+05          0.255154172          0.005804222
## Lag 2e+05          0.098559213         -0.004915787
## Lag 3e+05          0.037375858          0.002325456
## Lag 4e+05          0.017349115          0.015524794
## Lag 5e+05          0.005757827          0.022083184
## Chain 4 
##                   edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.0000000000          1.0000000000          1.0000000000
## Lag 1e+05  0.0107561813         -0.0005184371          0.0270551510
## Lag 2e+05 -0.0004782635         -0.0213638452          0.0009989092
## Lag 3e+05  0.0031789721         -0.0194695799          0.0170358783
## Lag 4e+05  0.0074754519         -0.0330261535          0.0198503988
## Lag 5e+05  0.0044960120          0.0115033137          0.0248903117
##           nodefactor.race..wa.B nodefactor.race..wa.H nodematch.race..wa.B
## Lag 0               1.000000000            1.00000000          1.000000000
## Lag 1e+05           0.030638882            0.04278515          0.072893477
## Lag 2e+05           0.018668703            0.03889139          0.008819404
## Lag 3e+05          -0.005133513            0.01089026          0.008439988
## Lag 4e+05          -0.004081987           -0.00390614          0.004266348
## Lag 5e+05           0.003129103           -0.01273688          0.016858548
##           nodematch.race..wa.H nodematch.race..wa.O
## Lag 0             1.0000000000          1.000000000
## Lag 1e+05         0.2375936227          0.007978804
## Lag 2e+05         0.0579818834         -0.003135887
## Lag 3e+05         0.0149737602         -0.004675465
## Lag 4e+05        -0.0004692792          0.007767258
## Lag 5e+05        -0.0164553055          0.005543742
## 
## Sample statistics burn-in diagnostic (Geweke):
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##                0.6597                0.8859                0.2569 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodematch.race..wa.B 
##               -0.4230               -1.8419               -0.5330 
##  nodematch.race..wa.H  nodematch.race..wa.O 
##               -1.5012                1.4016 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.50943011            0.37565741            0.79727829 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodematch.race..wa.B 
##            0.67230278            0.06549209            0.59405953 
##  nodematch.race..wa.H  nodematch.race..wa.O 
##            0.13330922            0.16104964 
## Joint P-value (lower = worse):  0.4008152 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##                1.4684               -0.4828                2.0007 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodematch.race..wa.B 
##                0.3025                1.5219                2.0254 
##  nodematch.race..wa.H  nodematch.race..wa.O 
##                1.4606                1.3057 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.14198484            0.62926000            0.04542308 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodematch.race..wa.B 
##            0.76226699            0.12802286            0.04282867 
##  nodematch.race..wa.H  nodematch.race..wa.O 
##            0.14412659            0.19165580 
## Joint P-value (lower = worse):  0.1298688 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##              -0.35143              -1.88365               1.06622 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodematch.race..wa.B 
##              -0.40138               1.02980              -0.09193 
##  nodematch.race..wa.H  nodematch.race..wa.O 
##               0.52247              -0.69290 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.72526297            0.05961224            0.28632547 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodematch.race..wa.B 
##            0.68814021            0.30310438            0.92675117 
##  nodematch.race..wa.H  nodematch.race..wa.O 
##            0.60134049            0.48837170 
## Joint P-value (lower = worse):  0.5479493 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##              -0.87791               0.09086               0.70552 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodematch.race..wa.B 
##              -1.20918              -0.86620               0.86729 
##  nodematch.race..wa.H  nodematch.race..wa.O 
##               1.38674               0.46245 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.3799926             0.9276013             0.4804885 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodematch.race..wa.B 
##             0.2265926             0.3863796             0.3857804 
##  nodematch.race..wa.H  nodematch.race..wa.O 
##             0.1655205             0.6437619 
## Joint P-value (lower = worse):  0.1683648 .
## Warning in formals(fun): argument is not a function

## 
## MCMC diagnostics shown here are from the last round of simulation, prior to computation of final parameter estimates. Because the final estimates are refinements of those used for this simulation run, these diagnostics may understate model performance. To directly assess the performance of the final model on in-model statistics, please use the GOF command: gof(ergmFitObject, GOF=~model).

Model 5

## Sample statistics summary:
## 
## Iterations = 1e+06:750900000
## Thinning interval = 1e+05 
## Number of chains = 4 
## Sample size per chain = 7500 
## 
## 1. Empirical mean and standard deviation for each variable,
##    plus standard error of the mean:
## 
##                           Mean     SD Naive SE Time-series SE
## edges                  0.30773 28.780  0.16616        0.16882
## nodefactor.deg.pers.1 -0.10860 17.687  0.10212        0.10212
## nodefactor.deg.pers.2  0.65457 18.613  0.10746        0.10706
## nodefactor.race..wa.B -0.02497 12.705  0.07335        0.07476
## nodefactor.race..wa.H -0.11593 17.477  0.10090        0.10925
## nodefactor.region.EW   0.25263 15.945  0.09206        0.09241
## nodefactor.region.OW   0.51080 29.500  0.17032        0.17027
## nodematch.race..wa.B  -0.07311  5.014  0.02895        0.03112
## nodematch.race..wa.H  -0.25268  8.743  0.05048        0.06592
## nodematch.race..wa.O   0.15164 26.418  0.15252        0.15432
## 
## 2. Quantiles for each variable:
## 
##                          2.5%     25%        50%    75%  97.5%
## edges                 -56.500 -19.500  5.000e-01 19.500 56.500
## nodefactor.deg.pers.1 -35.000 -12.000  5.684e-14 12.000 35.000
## nodefactor.deg.pers.2 -36.000 -12.000  0.000e+00 13.000 37.025
## nodefactor.race..wa.B -24.834  -8.834  1.664e-01  8.166 25.166
## nodefactor.race..wa.H -34.844 -11.844  1.560e-01 11.156 34.156
## nodefactor.region.EW  -30.561 -10.561  4.392e-01 11.439 31.439
## nodefactor.region.OW  -57.131 -19.131  8.694e-01 19.869 58.869
## nodematch.race..wa.B   -9.177  -3.177 -1.769e-01  2.823  9.823
## nodematch.race..wa.H  -17.300  -6.300 -3.003e-01  5.700 16.700
## nodematch.race..wa.O  -51.946 -17.946  5.397e-02 18.054 52.054
## 
## 
## Sample statistics cross-correlations:
##                           edges nodefactor.deg.pers.1
## edges                 1.0000000            0.40396471
## nodefactor.deg.pers.1 0.4039647            1.00000000
## nodefactor.deg.pers.2 0.4165221            0.04556771
## nodefactor.race..wa.B 0.2547034            0.11010520
## nodefactor.race..wa.H 0.3544556            0.16496778
## nodefactor.region.EW  0.3618646            0.14675518
## nodefactor.region.OW  0.6232773            0.23730057
## nodematch.race..wa.B  0.1060823            0.05009459
## nodematch.race..wa.H  0.1348891            0.06570727
## nodematch.race..wa.O  0.8188100            0.31943367
##                       nodefactor.deg.pers.2 nodefactor.race..wa.B
## edges                            0.41652210            0.25470343
## nodefactor.deg.pers.1            0.04556771            0.11010520
## nodefactor.deg.pers.2            1.00000000            0.11326212
## nodefactor.race..wa.B            0.11326212            1.00000000
## nodefactor.race..wa.H            0.15764390            0.04568274
## nodefactor.region.EW             0.16224203            0.04554084
## nodefactor.region.OW             0.25688763            0.12612004
## nodematch.race..wa.B             0.04423138            0.57492530
## nodematch.race..wa.H             0.05852404           -0.03842465
## nodematch.race..wa.O             0.33153153           -0.06020016
##                       nodefactor.race..wa.H nodefactor.region.EW
## edges                            0.35445562          0.361864612
## nodefactor.deg.pers.1            0.16496778          0.146755181
## nodefactor.deg.pers.2            0.15764390          0.162242030
## nodefactor.race..wa.B            0.04568274          0.045540841
## nodefactor.race..wa.H            1.00000000          0.237324037
## nodefactor.region.EW             0.23732404          1.000000000
## nodefactor.region.OW             0.20011369          0.052898411
## nodematch.race..wa.B            -0.01951699          0.004090968
## nodematch.race..wa.H             0.59195493          0.119488153
## nodematch.race..wa.O            -0.06956484          0.264134522
##                       nodefactor.region.OW nodematch.race..wa.B
## edges                           0.62327731          0.106082278
## nodefactor.deg.pers.1           0.23730057          0.050094593
## nodefactor.deg.pers.2           0.25688763          0.044231382
## nodefactor.race..wa.B           0.12612004          0.574925296
## nodefactor.race..wa.H           0.20011369         -0.019516987
## nodefactor.region.EW            0.05289841          0.004090968
## nodefactor.region.OW            1.00000000          0.049680598
## nodematch.race..wa.B            0.04968060          1.000000000
## nodematch.race..wa.H            0.07131174          0.022399286
## nodematch.race..wa.O            0.52993717          0.018386733
##                       nodematch.race..wa.H nodematch.race..wa.O
## edges                           0.13488906           0.81880999
## nodefactor.deg.pers.1           0.06570727           0.31943367
## nodefactor.deg.pers.2           0.05852404           0.33153153
## nodefactor.race..wa.B          -0.03842465          -0.06020016
## nodefactor.race..wa.H           0.59195493          -0.06956484
## nodefactor.region.EW            0.11948815           0.26413452
## nodefactor.region.OW            0.07131174           0.52993717
## nodematch.race..wa.B            0.02239929           0.01838673
## nodematch.race..wa.H            1.00000000           0.07491606
## nodematch.race..wa.O            0.07491606           1.00000000
## 
## Sample statistics auto-correlation:
## Chain 1 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05  0.008542287           0.002980698          -0.001923692
## Lag 2e+05 -0.009072005           0.009506289           0.001028801
## Lag 3e+05  0.014547358          -0.014283343           0.005353628
## Lag 4e+05 -0.007670068          -0.018298603           0.006005629
## Lag 5e+05 -0.004890381           0.009550699          -0.012004368
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0              1.0000000000          1.0000000000         1.0000000000
## Lag 1e+05          0.0208318802          0.0638244491         0.0123862384
## Lag 2e+05         -0.0128373565          0.0214285963        -0.0135528397
## Lag 3e+05          0.0115916808          0.0009343124         0.0110176195
## Lag 4e+05          0.0050129720         -0.0061646036        -0.0004921546
## Lag 5e+05         -0.0005475854         -0.0182401104        -0.0183672128
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000          1.000000000          1.000000000
## Lag 1e+05          0.018976014          0.097469770          0.254299214
## Lag 2e+05          0.022022062         -0.002958832          0.062433451
## Lag 3e+05          0.006651555         -0.005238543          0.007076696
## Lag 4e+05          0.002555040          0.005693267          0.011182032
## Lag 5e+05         -0.016968764          0.008590999          0.003805060
##           nodematch.race..wa.O
## Lag 0              1.000000000
## Lag 1e+05          0.009284854
## Lag 2e+05         -0.004895938
## Lag 3e+05          0.017894637
## Lag 4e+05          0.004252024
## Lag 5e+05         -0.003601692
## Chain 2 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05  0.012540947           0.013730624           0.030136490
## Lag 2e+05 -0.004804905          -0.001373605          -0.006692354
## Lag 3e+05 -0.006953238           0.001192709           0.004968828
## Lag 4e+05 -0.020140885          -0.008833628          -0.011103569
## Lag 5e+05  0.005872099          -0.003611982           0.003926455
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0              1.0000000000           1.000000000          1.000000000
## Lag 1e+05          0.0038987016           0.078017331          0.037281231
## Lag 2e+05         -0.0032975842           0.020904930          0.011249512
## Lag 3e+05          0.0025644764           0.005038043          0.003342170
## Lag 4e+05         -0.0057375707          -0.004827892         -0.011970714
## Lag 5e+05          0.0007974622           0.010332627         -0.006141835
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000          1.000000000          1.000000000
## Lag 1e+05          0.009044004          0.057341055          0.252366737
## Lag 2e+05          0.019061559         -0.003967935          0.079405384
## Lag 3e+05         -0.003010424          0.004075459          0.023944863
## Lag 4e+05         -0.019026916         -0.011249797          0.019379553
## Lag 5e+05         -0.002962861         -0.007594571          0.007173682
##           nodematch.race..wa.O
## Lag 0             1.0000000000
## Lag 1e+05         0.0101037678
## Lag 2e+05         0.0061156180
## Lag 3e+05         0.0067374804
## Lag 4e+05        -0.0198017091
## Lag 5e+05         0.0008365555
## Chain 3 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0     1.000000000           1.000000000           1.000000000
## Lag 1e+05 0.003862936          -0.007734902          -0.005954799
## Lag 2e+05 0.033001452           0.004657990          -0.003215551
## Lag 3e+05 0.010096485          -0.009778390          -0.040095934
## Lag 4e+05 0.002495206           0.008946910          -0.022964576
## Lag 5e+05 0.028192475           0.022115349           0.004208261
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0              1.0000000000           1.000000000          1.000000000
## Lag 1e+05          0.0109494472           0.063644821          0.014407612
## Lag 2e+05         -0.0142187949           0.024153762          0.005975160
## Lag 3e+05          0.0148181044           0.007620804         -0.010829447
## Lag 4e+05         -0.0006675172           0.013842050          0.003400889
## Lag 5e+05          0.0064483657          -0.001202179         -0.001259285
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0             1.0000000000          1.000000000           1.00000000
## Lag 1e+05         0.0002154924          0.074716651           0.27177817
## Lag 2e+05         0.0148608685         -0.011833517           0.08482142
## Lag 3e+05         0.0099464331          0.008370527           0.03873098
## Lag 4e+05         0.0082310620         -0.001350167           0.02513453
## Lag 5e+05         0.0127924170         -0.011160009           0.02036346
##           nodematch.race..wa.O
## Lag 0             1.0000000000
## Lag 1e+05         0.0001284225
## Lag 2e+05         0.0188360287
## Lag 3e+05         0.0122445986
## Lag 4e+05        -0.0049103393
## Lag 5e+05         0.0357362603
## Chain 4 
##                   edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.0000000000           1.000000000           1.000000000
## Lag 1e+05 -0.0135912469          -0.001684571           0.021632994
## Lag 2e+05  0.0105364432          -0.004941211           0.010307321
## Lag 3e+05 -0.0042666278          -0.003747846          -0.011002925
## Lag 4e+05 -0.0002665829          -0.007986214          -0.005327683
## Lag 5e+05  0.0055398314           0.015499904          -0.005590757
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000           1.000000000          1.000000000
## Lag 1e+05           0.033573454           0.062566637          0.008869042
## Lag 2e+05           0.021078872           0.025947791         -0.016056179
## Lag 3e+05          -0.008178511           0.004778550          0.008485616
## Lag 4e+05           0.020617062           0.012945109         -0.025130480
## Lag 5e+05          -0.007265961          -0.000330355         -0.001663378
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000          1.000000000          1.000000000
## Lag 1e+05         -0.017741248          0.074562374          0.248381670
## Lag 2e+05         -0.007331964          0.003051506          0.063443280
## Lag 3e+05          0.003760935         -0.005910881          0.024090969
## Lag 4e+05          0.007215346         -0.003921630          0.013958284
## Lag 5e+05         -0.030226059         -0.004940401         -0.002631113
##           nodematch.race..wa.O
## Lag 0             1.0000000000
## Lag 1e+05        -0.0197973178
## Lag 2e+05         0.0004861966
## Lag 3e+05        -0.0062239339
## Lag 4e+05         0.0189977818
## Lag 5e+05        -0.0065658652
## 
## Sample statistics burn-in diagnostic (Geweke):
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               1.08054               1.75554              -0.93014 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               1.07170               0.12092               0.35507 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##               0.81130               0.73966              -0.05306 
##  nodematch.race..wa.O 
##               0.71775 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.27990281            0.07916769            0.35229774 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##            0.28385460            0.90375165            0.72254060 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##            0.41719450            0.45950724            0.95768467 
##  nodematch.race..wa.O 
##            0.47291002 
## Joint P-value (lower = worse):  0.8201166 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               0.47719              -1.05450               0.54443 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##              -0.26835              -0.20247              -0.34744 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##               0.02862              -0.08727              -1.04427 
##  nodematch.race..wa.O 
##               0.12505 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.6332295             0.2916534             0.5861450 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##             0.7884303             0.8395525             0.7282581 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##             0.9771682             0.9304547             0.2963606 
##  nodematch.race..wa.O 
##             0.9004860 
## Joint P-value (lower = worse):  0.7460311 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               0.40033              -1.13108              -0.35552 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               1.57601               0.80942               0.01774 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##               1.45095               0.84036               0.47417 
##  nodematch.race..wa.O 
##              -0.63682 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.6889163             0.2580217             0.7222034 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##             0.1150227             0.4182748             0.9858476 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##             0.1467937             0.4007068             0.6353768 
##  nodematch.race..wa.O 
##             0.5242391 
## Joint P-value (lower = worse):  0.5058335 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##                0.7484                0.3772                1.0303 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##                1.4117                0.8763                0.6156 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##                0.8934                2.1229                1.1259 
##  nodematch.race..wa.O 
##               -0.1002 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.4542075             0.7060321             0.3028814 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##             0.1580386             0.3808809             0.5381495 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##             0.3716267             0.0337628             0.2602005 
##  nodematch.race..wa.O 
##             0.9201931 
## Joint P-value (lower = worse):  0.4846874 .
## Warning in formals(fun): argument is not a function

## 
## MCMC diagnostics shown here are from the last round of simulation, prior to computation of final parameter estimates. Because the final estimates are refinements of those used for this simulation run, these diagnostics may understate model performance. To directly assess the performance of the final model on in-model statistics, please use the GOF command: gof(ergmFitObject, GOF=~model).

Model 6

## Sample statistics summary:
## 
## Iterations = 1e+06:750900000
## Thinning interval = 1e+05 
## Number of chains = 4 
## Sample size per chain = 7500 
## 
## 1. Empirical mean and standard deviation for each variable,
##    plus standard error of the mean:
## 
##                          Mean     SD Naive SE Time-series SE
## edges                 -0.3362 28.655  0.16544        0.16702
## nodefactor.deg.pers.1 -0.4341 17.457  0.10079        0.10211
## nodefactor.deg.pers.2  0.1446 18.665  0.10776        0.11130
## nodefactor.race..wa.B -0.2209 12.626  0.07290        0.07902
## nodefactor.race..wa.H -0.0809 17.251  0.09960        0.12541
## nodefactor.region.EW   0.1399 16.051  0.09267        0.09872
## nodefactor.region.OW  -0.2902 29.410  0.16980        0.17710
## nodematch.race..wa.B  -0.1316  4.993  0.02882        0.03651
## nodematch.race..wa.H  -0.1758  8.651  0.04995        0.08025
## nodematch.race..wa.O  -0.2148 26.192  0.15122        0.15481
## absdiff.sqrt.age      -0.1908 28.553  0.16485        0.16688
## 
## 2. Quantiles for each variable:
## 
##                          2.5%     25%        50%    75%  97.5%
## edges                 -56.500 -19.500 -5.000e-01 19.500 56.500
## nodefactor.deg.pers.1 -35.000 -12.000  5.684e-14 11.000 34.000
## nodefactor.deg.pers.2 -36.000 -12.000  0.000e+00 13.000 37.000
## nodefactor.race..wa.B -24.834  -8.834 -8.336e-01  8.166 25.166
## nodefactor.race..wa.H -33.844 -11.844  1.560e-01 11.156 33.156
## nodefactor.region.EW  -31.561 -10.561  4.392e-01 11.439 31.439
## nodefactor.region.OW  -57.156 -20.131 -1.306e-01 19.869 57.869
## nodematch.race..wa.B   -9.177  -3.177 -1.769e-01  2.823  9.823
## nodematch.race..wa.H  -17.300  -6.300 -3.003e-01  5.700 16.700
## nodematch.race..wa.O  -50.946 -17.946  5.397e-02 17.054 51.054
## absdiff.sqrt.age      -55.102 -19.759 -2.802e-01 19.326 55.985
## 
## 
## Sample statistics cross-correlations:
##                           edges nodefactor.deg.pers.1
## edges                 1.0000000            0.40050192
## nodefactor.deg.pers.1 0.4005019            1.00000000
## nodefactor.deg.pers.2 0.4183869            0.04159989
## nodefactor.race..wa.B 0.2642528            0.10678321
## nodefactor.race..wa.H 0.3493018            0.15617993
## nodefactor.region.EW  0.3678966            0.14652313
## nodefactor.region.OW  0.6211976            0.23835336
## nodematch.race..wa.B  0.1082677            0.04757942
## nodematch.race..wa.H  0.1184630            0.05344195
## nodematch.race..wa.O  0.8219387            0.32337512
## absdiff.sqrt.age      0.5478571            0.22263633
##                       nodefactor.deg.pers.2 nodefactor.race..wa.B
## edges                            0.41838689            0.26425277
## nodefactor.deg.pers.1            0.04159989            0.10678321
## nodefactor.deg.pers.2            1.00000000            0.12405573
## nodefactor.race..wa.B            0.12405573            1.00000000
## nodefactor.race..wa.H            0.16635599            0.03824666
## nodefactor.region.EW             0.15807638            0.04720871
## nodefactor.region.OW             0.25274689            0.13277548
## nodematch.race..wa.B             0.05043611            0.56936634
## nodematch.race..wa.H             0.06154356           -0.05195565
## nodematch.race..wa.O             0.33159568           -0.04661283
## absdiff.sqrt.age                 0.21914700            0.14299754
##                       nodefactor.race..wa.H nodefactor.region.EW
## edges                            0.34930178           0.36789662
## nodefactor.deg.pers.1            0.15617993           0.14652313
## nodefactor.deg.pers.2            0.16635599           0.15807638
## nodefactor.race..wa.B            0.03824666           0.04720871
## nodefactor.race..wa.H            1.00000000           0.23798291
## nodefactor.region.EW             0.23798291           1.00000000
## nodefactor.region.OW             0.20222363           0.05508460
## nodematch.race..wa.B            -0.01793503           0.01239415
## nodematch.race..wa.H             0.59341279           0.12177764
## nodematch.race..wa.O            -0.06875286           0.27388681
## absdiff.sqrt.age                 0.19222431           0.20188113
##                       nodefactor.region.OW nodematch.race..wa.B
## edges                           0.62119761           0.10826773
## nodefactor.deg.pers.1           0.23835336           0.04757942
## nodefactor.deg.pers.2           0.25274689           0.05043611
## nodefactor.race..wa.B           0.13277548           0.56936634
## nodefactor.race..wa.H           0.20222363          -0.01793503
## nodefactor.region.EW            0.05508460           0.01239415
## nodefactor.region.OW            1.00000000           0.05226688
## nodematch.race..wa.B            0.05226688           1.00000000
## nodematch.race..wa.H            0.06943719           0.01124761
## nodematch.race..wa.O            0.52712791           0.01935570
## absdiff.sqrt.age                0.34312596           0.05798407
##                       nodematch.race..wa.H nodematch.race..wa.O
## edges                           0.11846300           0.82193870
## nodefactor.deg.pers.1           0.05344195           0.32337512
## nodefactor.deg.pers.2           0.06154356           0.33159568
## nodefactor.race..wa.B          -0.05195565          -0.04661283
## nodefactor.race..wa.H           0.59341279          -0.06875286
## nodefactor.region.EW            0.12177764           0.27388681
## nodefactor.region.OW            0.06943719           0.52712791
## nodematch.race..wa.B            0.01124761           0.01935570
## nodematch.race..wa.H            1.00000000           0.05971564
## nodematch.race..wa.O            0.05971564           1.00000000
## absdiff.sqrt.age                0.06557502           0.45041850
##                       absdiff.sqrt.age
## edges                       0.54785710
## nodefactor.deg.pers.1       0.22263633
## nodefactor.deg.pers.2       0.21914700
## nodefactor.race..wa.B       0.14299754
## nodefactor.race..wa.H       0.19222431
## nodefactor.region.EW        0.20188113
## nodefactor.region.OW        0.34312596
## nodematch.race..wa.B        0.05798407
## nodematch.race..wa.H        0.06557502
## nodematch.race..wa.O        0.45041850
## absdiff.sqrt.age            1.00000000
## 
## Sample statistics auto-correlation:
## Chain 1 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05  0.039202698           0.007405939           0.023281406
## Lag 2e+05  0.008849025          -0.016241057          -0.007577300
## Lag 3e+05  0.009454120           0.011486136          -0.002168793
## Lag 4e+05 -0.000657467          -0.001989726           0.028260015
## Lag 5e+05  0.002190588          -0.006470060           0.003226308
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000            1.00000000         1.0000000000
## Lag 1e+05           0.079319973            0.11608748         0.0318022430
## Lag 2e+05           0.017679860            0.06564404         0.0363968005
## Lag 3e+05           0.020526775            0.03970924        -0.0024831145
## Lag 4e+05          -0.006065356            0.02360628         0.0008340641
## Lag 5e+05          -0.004804690            0.02912619         0.0086500716
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000         1.0000000000           1.00000000
## Lag 1e+05          0.017453746         0.2209212668           0.38040379
## Lag 2e+05          0.008128531         0.0617716086           0.19338139
## Lag 3e+05          0.025604259         0.0114050369           0.10119218
## Lag 4e+05         -0.013240204         0.0007844741           0.06659826
## Lag 5e+05         -0.020024262        -0.0057973106           0.04366684
##           nodematch.race..wa.O absdiff.sqrt.age
## Lag 0              1.000000000      1.000000000
## Lag 1e+05          0.018983881      0.026557556
## Lag 2e+05         -0.005272583      0.014193569
## Lag 3e+05          0.004027528      0.006734619
## Lag 4e+05          0.003253652      0.009283374
## Lag 5e+05         -0.007907678     -0.001544995
## Chain 2 
##                   edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.0000000000           1.000000000           1.000000000
## Lag 1e+05  0.0078898745           0.028432815           0.027609495
## Lag 2e+05  0.0090811246          -0.005411520           0.029650546
## Lag 3e+05 -0.0007572838          -0.010038814           0.011420085
## Lag 4e+05  0.0262206879           0.005338408          -0.009432917
## Lag 5e+05  0.0075471619           0.010833606          -0.021532926
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000            1.00000000          1.000000000
## Lag 1e+05           0.045893461            0.10976373          0.022179122
## Lag 2e+05          -0.010005613            0.06334136          0.021275726
## Lag 3e+05          -0.003425241            0.03131515          0.002323357
## Lag 4e+05          -0.004788523            0.03895106          0.004044513
## Lag 5e+05          -0.013860879            0.02538180         -0.008549236
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000          1.000000000           1.00000000
## Lag 1e+05          0.017347717          0.198277945           0.37852164
## Lag 2e+05          0.021237604          0.059849669           0.19832981
## Lag 3e+05          0.005099190          0.016794947           0.10168638
## Lag 4e+05         -0.002910693          0.007856374           0.06690046
## Lag 5e+05         -0.005287947          0.014912908           0.04195017
##           nodematch.race..wa.O absdiff.sqrt.age
## Lag 0              1.000000000      1.000000000
## Lag 1e+05          0.018759377      0.013886888
## Lag 2e+05          0.015500809      0.004446377
## Lag 3e+05         -0.002887867      0.001593400
## Lag 4e+05          0.001572251      0.013226695
## Lag 5e+05         -0.009081279     -0.003021454
## Chain 3 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05 -0.019725203           0.022900432           0.018809493
## Lag 2e+05 -0.007216561           0.010756420          -0.010814419
## Lag 3e+05  0.023719388           0.016496840           0.004405406
## Lag 4e+05  0.003679112           0.013443043           0.017629107
## Lag 5e+05 -0.011392314           0.008385525          -0.011035933
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0              1.0000000000            1.00000000          1.000000000
## Lag 1e+05          0.0750644258            0.13367725          0.048554505
## Lag 2e+05          0.0291280462            0.05316322          0.004705089
## Lag 3e+05          0.0069559287            0.04823712          0.009972577
## Lag 4e+05          0.0135765204            0.02074794          0.031866319
## Lag 5e+05         -0.0007771856            0.03188090         -0.008587300
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000           1.00000000           1.00000000
## Lag 1e+05          0.024425796           0.20222605           0.39837462
## Lag 2e+05          0.013746716           0.07333533           0.19537954
## Lag 3e+05          0.019487803           0.04153427           0.10507243
## Lag 4e+05         -0.005074453           0.01703641           0.07188453
## Lag 5e+05         -0.005113679          -0.01362445           0.04614958
##           nodematch.race..wa.O absdiff.sqrt.age
## Lag 0              1.000000000     1.0000000000
## Lag 1e+05         -0.003554326    -0.0006306819
## Lag 2e+05          0.008618389     0.0067487308
## Lag 3e+05          0.004473344    -0.0133906729
## Lag 4e+05         -0.006428943    -0.0050674595
## Lag 5e+05         -0.009561283    -0.0030457780
## Chain 4 
##                   edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.0000000000            1.00000000           1.000000000
## Lag 1e+05  0.0009841128           -0.01130665           0.014637155
## Lag 2e+05  0.0031935598           -0.01059705           0.010565625
## Lag 3e+05  0.0106141369            0.01963256           0.021235479
## Lag 4e+05 -0.0193698090           -0.01595087           0.008823975
## Lag 5e+05  0.0284708951            0.01006575           0.018021259
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0              1.0000000000            1.00000000          1.000000000
## Lag 1e+05          0.0684207264            0.10751559          0.023631187
## Lag 2e+05          0.0319234240            0.04458401          0.026131589
## Lag 3e+05         -0.0009753907            0.03102584         -0.006319167
## Lag 4e+05         -0.0177870064            0.03682934          0.010573238
## Lag 5e+05         -0.0068226408            0.01377024          0.003874323
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000          1.000000000           1.00000000
## Lag 1e+05          0.038106049          0.198210076           0.37514744
## Lag 2e+05          0.005911342          0.073822406           0.19022610
## Lag 3e+05          0.028541063          0.026117720           0.10473242
## Lag 4e+05          0.004222651         -0.004035867           0.05992105
## Lag 5e+05          0.016565473          0.006168649           0.03145967
##           nodematch.race..wa.O absdiff.sqrt.age
## Lag 0              1.000000000      1.000000000
## Lag 1e+05          0.004465693      0.021688135
## Lag 2e+05          0.001783250     -0.004907913
## Lag 3e+05          0.026146559     -0.003419824
## Lag 4e+05         -0.004399005     -0.007028599
## Lag 5e+05          0.026131602      0.011143464
## 
## Sample statistics burn-in diagnostic (Geweke):
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##                0.2090                0.3562               -0.8867 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##                1.0885               -0.4525               -0.1697 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##                0.4274                0.3379               -1.1139 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##               -0.5515                0.9429 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.8344286             0.7217158             0.3752216 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##             0.2763737             0.6508924             0.8652068 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##             0.6690743             0.7354059             0.2653077 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##             0.5813154             0.3457085 
## Joint P-value (lower = worse):  0.8507749 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               1.06739              -0.58502              -0.02166 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               0.25208               1.76780               0.72173 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##              -0.03252               0.28317               1.62998 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##               0.35058              -1.31267 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.28579799            0.55853706            0.98271872 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##            0.80098239            0.07709386            0.47046081 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##            0.97405537            0.77704349            0.10310661 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##            0.72590015            0.18929422 
## Joint P-value (lower = worse):  0.3949642 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               1.66087              -0.03154               0.07853 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               1.67385               0.16334               0.02989 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##               1.28603              -0.58260              -0.40245 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##               0.81519               1.05444 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.09674034            0.97484059            0.93740335 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##            0.09415921            0.87025360            0.97615249 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##            0.19843303            0.56016035            0.68735224 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##            0.41496159            0.29168302 
## Joint P-value (lower = worse):  0.4327105 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               0.11764               0.53568              -0.54115 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               1.27794              -0.48613              -0.13664 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##              -0.49623               1.32883              -0.97710 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##               0.08461              -1.31020 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.9063553             0.5921799             0.5884014 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##             0.2012720             0.6268758             0.8913127 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##             0.6197329             0.1839041             0.3285193 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##             0.9325682             0.1901273 
## Joint P-value (lower = worse):  0.5245738 .
## Warning in formals(fun): argument is not a function

## 
## MCMC diagnostics shown here are from the last round of simulation, prior to computation of final parameter estimates. Because the final estimates are refinements of those used for this simulation run, these diagnostics may understate model performance. To directly assess the performance of the final model on in-model statistics, please use the GOF command: gof(ergmFitObject, GOF=~model).

Model 7

## Sample statistics summary:
## 
## Iterations = 1e+06:750900000
## Thinning interval = 1e+05 
## Number of chains = 4 
## Sample size per chain = 7500 
## 
## 1. Empirical mean and standard deviation for each variable,
##    plus standard error of the mean:
## 
##                          Mean     SD Naive SE Time-series SE
## edges                  0.3409 28.954  0.16716        0.20509
## nodefactor.deg.pers.1  0.4337 17.596  0.10159        0.11444
## nodefactor.deg.pers.2  0.0205 18.787  0.10847        0.13030
## nodefactor.race..wa.B -0.1459 12.554  0.07248        0.08937
## nodefactor.race..wa.H  0.8537 17.322  0.10001        0.17970
## nodefactor.region.EW   0.2217 18.288  0.10558        0.16969
## nodefactor.region.OW   0.6510 33.296  0.19223        0.20469
## nodematch.race..wa.B  -0.1799  4.931  0.02847        0.04520
## nodematch.race..wa.H   0.1733  8.610  0.04971        0.13683
## nodematch.race..wa.O  -0.3078 26.309  0.15190        0.17291
## absdiff.sqrt.age      -0.1465 28.569  0.16494        0.17867
## 
## 2. Quantiles for each variable:
## 
##                          2.5%     25%        50%    75%  97.5%
## edges                 -55.500 -19.500  5.000e-01 19.500 57.500
## nodefactor.deg.pers.1 -34.000 -11.000  5.684e-14 12.000 35.000
## nodefactor.deg.pers.2 -37.000 -13.000  0.000e+00 13.000 37.000
## nodefactor.race..wa.B -24.834  -8.834  1.664e-01  8.166 24.191
## nodefactor.race..wa.H -32.844 -10.844  1.156e+00 12.156 35.156
## nodefactor.region.EW  -35.561 -11.561  4.392e-01 12.439 36.439
## nodefactor.region.OW  -64.131 -22.131 -1.306e-01 23.869 65.869
## nodematch.race..wa.B   -9.177  -3.177 -1.769e-01  2.823  9.823
## nodematch.race..wa.H  -16.300  -5.300 -3.003e-01  5.700 16.700
## nodematch.race..wa.O  -51.946 -17.946 -9.460e-01 17.054 51.054
## absdiff.sqrt.age      -55.139 -19.522 -4.877e-01 18.963 56.752
## 
## 
## Sample statistics cross-correlations:
##                           edges nodefactor.deg.pers.1
## edges                 1.0000000            0.39587037
## nodefactor.deg.pers.1 0.3958704            1.00000000
## nodefactor.deg.pers.2 0.4324802            0.04604631
## nodefactor.race..wa.B 0.2699487            0.10578419
## nodefactor.race..wa.H 0.3548353            0.16164935
## nodefactor.region.EW  0.3161591            0.13435631
## nodefactor.region.OW  0.5747786            0.20969378
## nodematch.race..wa.B  0.1041225            0.04363221
## nodematch.race..wa.H  0.1331056            0.07335396
## nodematch.race..wa.O  0.8228300            0.31832444
## absdiff.sqrt.age      0.5517659            0.21343431
##                       nodefactor.deg.pers.2 nodefactor.race..wa.B
## edges                            0.43248025            0.26994868
## nodefactor.deg.pers.1            0.04604631            0.10578419
## nodefactor.deg.pers.2            1.00000000            0.13289750
## nodefactor.race..wa.B            0.13289750            1.00000000
## nodefactor.race..wa.H            0.15813675            0.04815951
## nodefactor.region.EW             0.14796254            0.03673300
## nodefactor.region.OW             0.24118703            0.11877553
## nodematch.race..wa.B             0.05976305            0.56013482
## nodematch.race..wa.H             0.05531873           -0.04273603
## nodematch.race..wa.O             0.34736867           -0.04290010
## absdiff.sqrt.age                 0.23686395            0.15504950
##                       nodefactor.race..wa.H nodefactor.region.EW
## edges                            0.35483532          0.316159143
## nodefactor.deg.pers.1            0.16164935          0.134356313
## nodefactor.deg.pers.2            0.15813675          0.147962541
## nodefactor.race..wa.B            0.04815951          0.036732997
## nodefactor.race..wa.H            1.00000000          0.224770402
## nodefactor.region.EW             0.22477040          1.000000000
## nodefactor.region.OW             0.19096190         -0.006425365
## nodematch.race..wa.B            -0.01382299          0.013623290
## nodematch.race..wa.H             0.59690387          0.115037527
## nodematch.race..wa.O            -0.06482711          0.227290315
## absdiff.sqrt.age                 0.19812770          0.169504301
##                       nodefactor.region.OW nodematch.race..wa.B
## edges                          0.574778596          0.104122539
## nodefactor.deg.pers.1          0.209693779          0.043632214
## nodefactor.deg.pers.2          0.241187035          0.059763046
## nodefactor.race..wa.B          0.118775526          0.560134817
## nodefactor.race..wa.H          0.190961902         -0.013822993
## nodefactor.region.EW          -0.006425365          0.013623290
## nodefactor.region.OW           1.000000000          0.033492389
## nodematch.race..wa.B           0.033492389          1.000000000
## nodematch.race..wa.H           0.069243414          0.007052211
## nodematch.race..wa.O           0.489451963          0.014272111
## absdiff.sqrt.age               0.319325303          0.056101096
##                       nodematch.race..wa.H nodematch.race..wa.O
## edges                          0.133105552           0.82283002
## nodefactor.deg.pers.1          0.073353958           0.31832444
## nodefactor.deg.pers.2          0.055318733           0.34736867
## nodefactor.race..wa.B         -0.042736032          -0.04290010
## nodefactor.race..wa.H          0.596903872          -0.06482711
## nodefactor.region.EW           0.115037527           0.22729032
## nodefactor.region.OW           0.069243414           0.48945196
## nodematch.race..wa.B           0.007052211           0.01427211
## nodematch.race..wa.H           1.000000000           0.06765965
## nodematch.race..wa.O           0.067659652           1.00000000
## absdiff.sqrt.age               0.073090025           0.45050529
##                       absdiff.sqrt.age
## edges                       0.55176595
## nodefactor.deg.pers.1       0.21343431
## nodefactor.deg.pers.2       0.23686395
## nodefactor.race..wa.B       0.15504950
## nodefactor.race..wa.H       0.19812770
## nodefactor.region.EW        0.16950430
## nodefactor.region.OW        0.31932530
## nodematch.race..wa.B        0.05610110
## nodematch.race..wa.H        0.07309002
## nodematch.race..wa.O        0.45050529
## absdiff.sqrt.age            1.00000000
## 
## Sample statistics auto-correlation:
## Chain 1 
##                edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0     1.00000000           1.000000000            1.00000000
## Lag 1e+05 0.08826008           0.103947566            0.10704295
## Lag 2e+05 0.03503032           0.018560359            0.03972203
## Lag 3e+05 0.03054123           0.009470079            0.04943772
## Lag 4e+05 0.02189859           0.010352221            0.02718131
## Lag 5e+05 0.01461508           0.029428779            0.01596672
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000            1.00000000           1.00000000
## Lag 1e+05           0.118028260            0.28786662           0.32837545
## Lag 2e+05           0.050335796            0.16881826           0.17663378
## Lag 3e+05           0.037921167            0.13253545           0.10999496
## Lag 4e+05           0.041816160            0.08910375           0.06383326
## Lag 5e+05           0.007116593            0.05993399           0.03149067
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000           1.00000000            1.0000000
## Lag 1e+05          0.048887292           0.39016659            0.6749795
## Lag 2e+05          0.010883240           0.18136615            0.5031622
## Lag 3e+05          0.001000355           0.11887792            0.3914737
## Lag 4e+05          0.004674848           0.08217766            0.3016400
## Lag 5e+05          0.014585707           0.04875591            0.2443631
##           nodematch.race..wa.O absdiff.sqrt.age
## Lag 0              1.000000000     1.0000000000
## Lag 1e+05          0.088430011     0.0758183312
## Lag 2e+05          0.045131710     0.0336485313
## Lag 3e+05          0.022468001     0.0148307901
## Lag 4e+05          0.021031381     0.0009786120
## Lag 5e+05          0.009106743    -0.0005854615
## Chain 2 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000            1.00000000            1.00000000
## Lag 1e+05  0.078334705            0.09654722            0.11924434
## Lag 2e+05  0.043876068            0.03911747            0.05863247
## Lag 3e+05  0.029367331            0.02512674            0.03349835
## Lag 4e+05 -0.001398964            0.01349570            0.02964573
## Lag 5e+05  0.004757820            0.02119597            0.02799589
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000            1.00000000           1.00000000
## Lag 1e+05           0.134166576            0.29499482           0.32474530
## Lag 2e+05           0.050257453            0.17924878           0.18537783
## Lag 3e+05           0.024493139            0.14105384           0.10899688
## Lag 4e+05           0.009509095            0.08968222           0.06078367
## Lag 5e+05           0.018946917            0.07475333           0.06045068
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000           1.00000000            1.0000000
## Lag 1e+05          0.047777044           0.34932355            0.6657985
## Lag 2e+05          0.015904851           0.17192681            0.4914827
## Lag 3e+05          0.002498752           0.09539789            0.3721520
## Lag 4e+05          0.008040451           0.04346165            0.2877381
## Lag 5e+05         -0.006524582           0.02081922            0.2389953
##           nodematch.race..wa.O absdiff.sqrt.age
## Lag 0              1.000000000      1.000000000
## Lag 1e+05          0.086828290      0.055998143
## Lag 2e+05          0.059812064      0.014921365
## Lag 3e+05          0.027213278      0.009769602
## Lag 4e+05          0.004906727     -0.014839070
## Lag 5e+05          0.010804301     -0.005497228
## Chain 3 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000          1.0000000000
## Lag 1e+05  0.068598157           0.090574353          0.1038017951
## Lag 2e+05  0.035277835           0.050342182          0.0462273400
## Lag 3e+05  0.023009524           0.019999030          0.0232689254
## Lag 4e+05  0.009775106           0.016609135          0.0061996542
## Lag 5e+05 -0.008529641           0.001941454          0.0002037566
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000            1.00000000           1.00000000
## Lag 1e+05           0.147232131            0.27282434           0.31171121
## Lag 2e+05           0.047006967            0.17755842           0.15860930
## Lag 3e+05           0.029249513            0.11873364           0.08303764
## Lag 4e+05           0.031008926            0.09196339           0.06792686
## Lag 5e+05           0.005782397            0.08422395           0.05603013
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0             1.0000000000           1.00000000            1.0000000
## Lag 1e+05         0.0748416581           0.35939108            0.6628109
## Lag 2e+05         0.0038072322           0.16013384            0.4923825
## Lag 3e+05        -0.0016988172           0.09687858            0.3846767
## Lag 4e+05         0.0005295012           0.05966376            0.3098145
## Lag 5e+05         0.0121711277           0.02564017            0.2543966
##           nodematch.race..wa.O absdiff.sqrt.age
## Lag 0              1.000000000      1.000000000
## Lag 1e+05          0.075998145      0.062445629
## Lag 2e+05          0.038242493      0.019762294
## Lag 3e+05          0.022165940      0.007175586
## Lag 4e+05          0.008042402     -0.003784013
## Lag 5e+05         -0.011723425      0.001046195
## Chain 4 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0     1.0000000000          1.0000000000            1.00000000
## Lag 1e+05 0.0716115758          0.0904722936            0.12801040
## Lag 2e+05 0.0293412587          0.0110793718            0.04990259
## Lag 3e+05 0.0205640434          0.0097471058            0.01885096
## Lag 4e+05 0.0204827130          0.0220562275            0.01091646
## Lag 5e+05 0.0006692579          0.0004111998            0.01335019
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000            1.00000000           1.00000000
## Lag 1e+05           0.161623731            0.28408805           0.30930781
## Lag 2e+05           0.094821239            0.17671179           0.15268955
## Lag 3e+05           0.035337216            0.12588104           0.09408452
## Lag 4e+05           0.018497341            0.10969674           0.06068108
## Lag 5e+05          -0.009187217            0.09590119           0.05227518
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000          1.000000000            1.0000000
## Lag 1e+05          0.059269793          0.373886675            0.6615745
## Lag 2e+05          0.022089126          0.183853153            0.4894539
## Lag 3e+05         -0.007940382          0.097552211            0.3824800
## Lag 4e+05          0.005960086          0.041712722            0.3127496
## Lag 5e+05         -0.001017895          0.002672131            0.2546519
##           nodematch.race..wa.O absdiff.sqrt.age
## Lag 0              1.000000000      1.000000000
## Lag 1e+05          0.082007134      0.058243964
## Lag 2e+05          0.030392459      0.039378110
## Lag 3e+05          0.009830594      0.011875120
## Lag 4e+05          0.015019353      0.018135648
## Lag 5e+05          0.003382568     -0.002697232
## 
## Sample statistics burn-in diagnostic (Geweke):
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##                1.3681                0.9388                0.1340 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##                0.4567                0.4803                1.4241 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##                0.3031                0.8692                0.5920 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##                1.3869                1.6797 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.17127429            0.34783633            0.89340469 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##            0.64789914            0.63103004            0.15442436 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##            0.76181953            0.38476380            0.55386245 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##            0.16548182            0.09301537 
## Joint P-value (lower = worse):  0.8846388 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##              -0.27343               0.47075               0.26208 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               0.27454              -0.34724               0.79363 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##              -0.14319              -0.09275              -1.00121 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##              -1.00861              -0.47683 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.7845229             0.6378182             0.7932576 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##             0.7836676             0.7284113             0.4274134 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##             0.8861427             0.9261030             0.3167276 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##             0.3131636             0.6334858 
## Joint P-value (lower = worse):  0.3837697 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               0.72014               0.91675               0.37306 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##              -0.72357              -0.09954              -0.66263 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##               1.36336              -1.25520              -1.27408 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##               0.48043               0.34133 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.4714368             0.3592718             0.7091049 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##             0.4693275             0.9207113             0.5075692 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##             0.1727696             0.2094058             0.2026350 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##             0.6309237             0.7328565 
## Joint P-value (lower = worse):  0.4357829 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##              -0.55136               0.59264               0.64923 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##              -0.07306              -1.24957               0.61222 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##              -0.68118               1.22486               0.22372 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##               0.51968              -1.59813 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##             0.5813838             0.5534193             0.5161887 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##             0.9417601             0.2114577             0.5403951 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##             0.4957571             0.2206287             0.8229732 
##  nodematch.race..wa.O      absdiff.sqrt.age 
##             0.6032836             0.1100151 
## Joint P-value (lower = worse):  0.3235072 .
## Warning in formals(fun): argument is not a function

## 
## MCMC diagnostics shown here are from the last round of simulation, prior to computation of final parameter estimates. Because the final estimates are refinements of those used for this simulation run, these diagnostics may understate model performance. To directly assess the performance of the final model on in-model statistics, please use the GOF command: gof(ergmFitObject, GOF=~model).

Model 8

## Sample statistics summary:
## 
## Iterations = 1e+06:750900000
## Thinning interval = 1e+05 
## Number of chains = 4 
## Sample size per chain = 7500 
## 
## 1. Empirical mean and standard deviation for each variable,
##    plus standard error of the mean:
## 
##                          Mean     SD Naive SE Time-series SE
## edges                  0.6674 28.719  0.16581         0.1819
## nodefactor.deg.pers.1  0.1559 17.576  0.10147         0.1109
## nodefactor.deg.pers.2 -1.1388 18.733  0.10816         0.1225
## nodefactor.race..wa.B  0.2557 12.532  0.07235         0.0848
## nodefactor.race..wa.H -0.6082 17.367  0.10027         0.1651
## nodefactor.region.EW   0.1618 17.769  0.10259         0.1490
## nodefactor.region.OW  -0.2508 32.507  0.18768         0.2036
## nodematch.race..wa.B   0.1578  4.988  0.02880         0.0429
## nodematch.race..wa.H  -0.4796  8.768  0.05062         0.1212
## nodematch.race..wa.O   0.8242 26.327  0.15200         0.1685
## absdiff.sqrt.age       0.4665 28.286  0.16331         0.1742
## nodematch.region       0.1894 29.297  0.16915         0.1904
## 
## 2. Quantiles for each variable:
## 
##                          2.5%     25%        50%    75%  97.5%
## edges                 -55.500 -18.500  5.000e-01 20.500 56.500
## nodefactor.deg.pers.1 -34.000 -12.000  5.684e-14 12.000 35.000
## nodefactor.deg.pers.2 -37.000 -14.000 -1.000e+00 11.000 35.000
## nodefactor.race..wa.B -23.834  -7.834  1.664e-01  9.166 25.166
## nodefactor.race..wa.H -34.844 -11.844 -8.440e-01 11.156 33.156
## nodefactor.region.EW  -34.561 -11.561  4.392e-01 12.439 34.439
## nodefactor.region.OW  -64.131 -22.131 -1.306e-01 21.869 63.869
## nodematch.race..wa.B   -9.177  -3.177 -1.769e-01  3.823  9.823
## nodematch.race..wa.H  -17.300  -6.300 -3.003e-01  5.700 16.700
## nodematch.race..wa.O  -50.946 -16.946  1.054e+00 19.054 52.054
## absdiff.sqrt.age      -54.693 -18.678  3.171e-01 19.347 56.226
## nodematch.region      -57.450 -19.450 -4.500e-01 19.550 57.550
## 
## 
## Sample statistics cross-correlations:
##                           edges nodefactor.deg.pers.1
## edges                 1.0000000            0.39882373
## nodefactor.deg.pers.1 0.3988237            1.00000000
## nodefactor.deg.pers.2 0.4253927            0.04655527
## nodefactor.race..wa.B 0.2543469            0.11224116
## nodefactor.race..wa.H 0.3522774            0.15028060
## nodefactor.region.EW  0.3272193            0.12794292
## nodefactor.region.OW  0.5760714            0.21740494
## nodematch.race..wa.B  0.1057674            0.05401513
## nodematch.race..wa.H  0.1276504            0.05001325
## nodematch.race..wa.O  0.8222503            0.31871374
## absdiff.sqrt.age      0.5447800            0.20999808
## nodematch.region      0.8767400            0.35326691
##                       nodefactor.deg.pers.2 nodefactor.race..wa.B
## edges                            0.42539265            0.25434693
## nodefactor.deg.pers.1            0.04655527            0.11224116
## nodefactor.deg.pers.2            1.00000000            0.12553490
## nodefactor.race..wa.B            0.12553490            1.00000000
## nodefactor.race..wa.H            0.16807656            0.04419428
## nodefactor.region.EW             0.15086021            0.02440577
## nodefactor.region.OW             0.22873429            0.11507395
## nodematch.race..wa.B             0.05478773            0.57149448
## nodematch.race..wa.H             0.07324067           -0.03950840
## nodematch.race..wa.O             0.33995698           -0.05523210
## absdiff.sqrt.age                 0.23306318            0.13995024
## nodematch.region                 0.36982364            0.22943954
##                       nodefactor.race..wa.H nodefactor.region.EW
## edges                           0.352277377           0.32721929
## nodefactor.deg.pers.1           0.150280602           0.12794292
## nodefactor.deg.pers.2           0.168076555           0.15086021
## nodefactor.race..wa.B           0.044194283           0.02440577
## nodefactor.race..wa.H           1.000000000           0.24556499
## nodefactor.region.EW            0.245564986           1.00000000
## nodefactor.region.OW            0.189217659           0.01256993
## nodematch.race..wa.B           -0.003062516           0.01471581
## nodematch.race..wa.H            0.598192290           0.13931629
## nodematch.race..wa.O           -0.065845028           0.23580024
## absdiff.sqrt.age                0.194141292           0.16795694
## nodematch.region                0.303595438           0.23686343
##                       nodefactor.region.OW nodematch.race..wa.B
## edges                           0.57607137          0.105767360
## nodefactor.deg.pers.1           0.21740494          0.054015129
## nodefactor.deg.pers.2           0.22873429          0.054787727
## nodefactor.race..wa.B           0.11507395          0.571494482
## nodefactor.race..wa.H           0.18921766         -0.003062516
## nodefactor.region.EW            0.01256993          0.014715813
## nodefactor.region.OW            1.00000000          0.034407436
## nodematch.race..wa.B            0.03440744          1.000000000
## nodematch.race..wa.H            0.06539252          0.032055927
## nodematch.race..wa.O            0.48913232          0.015016380
## absdiff.sqrt.age                0.31245773          0.052566292
## nodematch.region                0.49090371          0.102124543
##                       nodematch.race..wa.H nodematch.race..wa.O
## edges                           0.12765039           0.82225030
## nodefactor.deg.pers.1           0.05001325           0.31871374
## nodefactor.deg.pers.2           0.07324067           0.33995698
## nodefactor.race..wa.B          -0.03950840          -0.05523210
## nodefactor.race..wa.H           0.59819229          -0.06584503
## nodefactor.region.EW            0.13931629           0.23580024
## nodefactor.region.OW            0.06539252           0.48913232
## nodematch.race..wa.B            0.03205593           0.01501638
## nodematch.race..wa.H            1.00000000           0.06833750
## nodematch.race..wa.O            0.06833750           1.00000000
## absdiff.sqrt.age                0.06123585           0.44211889
## nodematch.region                0.11351031           0.72161597
##                       absdiff.sqrt.age nodematch.region
## edges                       0.54478003        0.8767400
## nodefactor.deg.pers.1       0.20999808        0.3532669
## nodefactor.deg.pers.2       0.23306318        0.3698236
## nodefactor.race..wa.B       0.13995024        0.2294395
## nodefactor.race..wa.H       0.19414129        0.3035954
## nodefactor.region.EW        0.16795694        0.2368634
## nodefactor.region.OW        0.31245773        0.4909037
## nodematch.race..wa.B        0.05256629        0.1021245
## nodematch.race..wa.H        0.06123585        0.1135103
## nodematch.race..wa.O        0.44211889        0.7216160
## absdiff.sqrt.age            1.00000000        0.4790433
## nodematch.region            0.47904330        1.0000000
## 
## Sample statistics auto-correlation:
## Chain 1 
##                edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0     1.00000000           1.000000000            1.00000000
## Lag 1e+05 0.06319881           0.042584481            0.07785070
## Lag 2e+05 0.02901302           0.047480914            0.02300900
## Lag 3e+05 0.02957396           0.017559286            0.01665821
## Lag 4e+05 0.03339324          -0.004741893            0.02759303
## Lag 5e+05 0.02172360          -0.003271898            0.01812706
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0              1.0000000000            1.00000000           1.00000000
## Lag 1e+05          0.1118295831            0.20460716           0.22671890
## Lag 2e+05          0.0468231260            0.12723742           0.11488121
## Lag 3e+05          0.0130329440            0.09454750           0.06397098
## Lag 4e+05          0.0109854321            0.08078360           0.04484796
## Lag 5e+05         -0.0008268687            0.06205369           0.02046217
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000          1.000000000            1.0000000
## Lag 1e+05          0.070771600          0.329252257            0.5505790
## Lag 2e+05          0.009957765          0.149025650            0.3970296
## Lag 3e+05          0.024655149          0.070483496            0.2910314
## Lag 4e+05          0.019105841          0.024784916            0.2226306
## Lag 5e+05          0.021343673          0.009288153            0.1640154
##           nodematch.race..wa.O absdiff.sqrt.age nodematch.region
## Lag 0               1.00000000      1.000000000      1.000000000
## Lag 1e+05           0.08806448      0.058898532      0.095614874
## Lag 2e+05           0.02941725      0.019374321      0.047164713
## Lag 3e+05           0.02401473      0.018724482      0.045279657
## Lag 4e+05           0.01939618      0.009613618      0.029632498
## Lag 5e+05           0.01001271     -0.002262934      0.009240833
## Chain 2 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05  0.052379418           0.076003822           0.076860368
## Lag 2e+05  0.005765307           0.039212090           0.032646476
## Lag 3e+05  0.002331449           0.008994256           0.022003325
## Lag 4e+05  0.006942833           0.018363512           0.001016604
## Lag 5e+05 -0.008699272           0.007619185           0.002431188
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000            1.00000000           1.00000000
## Lag 1e+05           0.103820773            0.22215297           0.23049924
## Lag 2e+05           0.039564130            0.14324235           0.10583693
## Lag 3e+05           0.025353703            0.08993701           0.04268079
## Lag 4e+05           0.007950401            0.06286781           0.02462168
## Lag 5e+05          -0.002109837            0.05969175           0.02807392
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000           1.00000000            1.0000000
## Lag 1e+05          0.051972525           0.33115994            0.5506864
## Lag 2e+05         -0.014782457           0.14637483            0.3857423
## Lag 3e+05          0.009754428           0.09126691            0.2757427
## Lag 4e+05          0.011911817           0.05270805            0.2205706
## Lag 5e+05         -0.007504480           0.02147032            0.1892104
##           nodematch.race..wa.O absdiff.sqrt.age nodematch.region
## Lag 0              1.000000000      1.000000000      1.000000000
## Lag 1e+05          0.062345937      0.055122278      0.080214928
## Lag 2e+05          0.008739686      0.020940225      0.027256434
## Lag 3e+05         -0.001369010      0.011490805      0.003234749
## Lag 4e+05          0.014468063     -0.006033567      0.004398724
## Lag 5e+05         -0.018166102      0.006446410     -0.013199417
## Chain 3 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0     1.000000000           1.000000000           1.000000000
## Lag 1e+05 0.050490994           0.044537263           0.084113753
## Lag 2e+05 0.018789794           0.010959563           0.033040972
## Lag 3e+05 0.005876605           0.001727206           0.016975710
## Lag 4e+05 0.010477310           0.024201449           0.009136411
## Lag 5e+05 0.015358555           0.018882094           0.008408206
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0                1.00000000            1.00000000           1.00000000
## Lag 1e+05            0.11557796            0.21432242           0.24115464
## Lag 2e+05            0.05130032            0.15373107           0.12612823
## Lag 3e+05            0.02744415            0.11271971           0.09335819
## Lag 4e+05            0.01548805            0.09403663           0.04582351
## Lag 5e+05            0.01140012            0.08150393           0.03813721
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0              1.000000000           1.00000000            1.0000000
## Lag 1e+05          0.040277728           0.33768176            0.5381494
## Lag 2e+05          0.010837205           0.15275135            0.3842432
## Lag 3e+05          0.008247860           0.06316776            0.3040152
## Lag 4e+05         -0.001113046           0.04116602            0.2379443
## Lag 5e+05          0.028986054           0.01940220            0.1876097
##           nodematch.race..wa.O absdiff.sqrt.age nodematch.region
## Lag 0             1.0000000000      1.000000000      1.000000000
## Lag 1e+05         0.0617929383      0.054530046      0.066595706
## Lag 2e+05         0.0074123385      0.016204495      0.031566977
## Lag 3e+05         0.0006515408      0.002208614      0.015593173
## Lag 4e+05         0.0045486721      0.005560654      0.005262058
## Lag 5e+05         0.0109524308      0.006042239      0.023499797
## Chain 4 
##                  edges nodefactor.deg.pers.1 nodefactor.deg.pers.2
## Lag 0      1.000000000           1.000000000           1.000000000
## Lag 1e+05  0.058938044           0.091086164           0.077650217
## Lag 2e+05  0.041301746           0.026739394           0.037998745
## Lag 3e+05 -0.005129468           0.009130540           0.027620863
## Lag 4e+05  0.007709218           0.007673822           0.001098832
## Lag 5e+05 -0.009760663           0.008534006           0.014686789
##           nodefactor.race..wa.B nodefactor.race..wa.H nodefactor.region.EW
## Lag 0               1.000000000            1.00000000           1.00000000
## Lag 1e+05           0.133514131            0.21555913           0.21495067
## Lag 2e+05           0.072984982            0.14557414           0.11696602
## Lag 3e+05           0.026534158            0.10804051           0.06277538
## Lag 4e+05           0.015418246            0.08753610           0.03396595
## Lag 5e+05          -0.009185847            0.08665772           0.03876947
##           nodefactor.region.OW nodematch.race..wa.B nodematch.race..wa.H
## Lag 0             1.0000000000           1.00000000            1.0000000
## Lag 1e+05         0.0463043805           0.33230131            0.5560083
## Lag 2e+05         0.0229144770           0.14936175            0.4112324
## Lag 3e+05        -0.0056222946           0.07214151            0.3096621
## Lag 4e+05        -0.0179520002           0.02874199            0.2377860
## Lag 5e+05        -0.0006059438           0.01928351            0.1924291
##           nodematch.race..wa.O absdiff.sqrt.age nodematch.region
## Lag 0              1.000000000      1.000000000      1.000000000
## Lag 1e+05          0.076224830      0.033193761      0.083881016
## Lag 2e+05          0.047461216      0.027590533      0.046255084
## Lag 3e+05          0.009418221      0.005059686      0.017409631
## Lag 4e+05          0.031811648      0.012856321      0.020960200
## Lag 5e+05          0.006959893     -0.005689170      0.006420059
## 
## Sample statistics burn-in diagnostic (Geweke):
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               -0.9124               -0.4414                0.3485 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##                0.6907               -0.6099                0.7500 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##               -2.0339                1.4853               -1.2363 
##  nodematch.race..wa.O      absdiff.sqrt.age      nodematch.region 
##               -1.6250                0.3673               -0.7370 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.36154374            0.65890712            0.72748106 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##            0.48973801            0.54195933            0.45326564 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##            0.04195738            0.13746165            0.21633351 
##  nodematch.race..wa.O      absdiff.sqrt.age      nodematch.region 
##            0.10417130            0.71342374            0.46114374 
## Joint P-value (lower = worse):  0.1147244 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##                1.1884                0.5781                0.8807 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               -0.2619                0.8986                2.2547 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##                1.0678               -0.7907                1.0749 
##  nodematch.race..wa.O      absdiff.sqrt.age      nodematch.region 
##                1.4659               -0.4074                0.8391 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.23465944            0.56320956            0.37849606 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##            0.79337070            0.36886153            0.02415261 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##            0.28562255            0.42909614            0.28242800 
##  nodematch.race..wa.O      absdiff.sqrt.age      nodematch.region 
##            0.14268622            0.68372391            0.40138684 
## Joint P-value (lower = worse):  0.4359832 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               -1.0561               -0.3869                0.2332 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               -1.0409               -0.6943                0.3325 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##               -1.6928               -0.5937               -0.5064 
##  nodematch.race..wa.O      absdiff.sqrt.age      nodematch.region 
##               -0.5123                0.4037                0.5283 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##            0.29091814            0.69884542            0.81562757 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##            0.29790232            0.48750681            0.73951277 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##            0.09048624            0.55271232            0.61260592 
##  nodematch.race..wa.O      absdiff.sqrt.age      nodematch.region 
##            0.60847016            0.68646954            0.59732331 
## Joint P-value (lower = worse):  0.2496046 .
## Chain 4 
## 
## Fraction in 1st window = 0.1
## Fraction in 2nd window = 0.5 
## 
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##               -1.0081               -0.7275               -1.0099 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##               -2.2469               -0.6343                0.3251 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##                1.3355               -2.7636               -1.3202 
##  nodematch.race..wa.O      absdiff.sqrt.age      nodematch.region 
##               -0.7380                0.9269               -0.8632 
## 
## Individual P-values (lower = worse):
##                 edges nodefactor.deg.pers.1 nodefactor.deg.pers.2 
##           0.313388155           0.466939175           0.312523915 
## nodefactor.race..wa.B nodefactor.race..wa.H  nodefactor.region.EW 
##           0.024649387           0.525875683           0.745113623 
##  nodefactor.region.OW  nodematch.race..wa.B  nodematch.race..wa.H 
##           0.181700411           0.005716529           0.186761457 
##  nodematch.race..wa.O      absdiff.sqrt.age      nodematch.region 
##           0.460528190           0.353956081           0.388052612 
## Joint P-value (lower = worse):  0.07631727 .
## Warning in formals(fun): argument is not a function

## 
## MCMC diagnostics shown here are from the last round of simulation, prior to computation of final parameter estimates. Because the final estimates are refinements of those used for this simulation run, these diagnostics may understate model performance. To directly assess the performance of the final model on in-model statistics, please use the GOF command: gof(ergmFitObject, GOF=~model).

Summary of model fit

Model 1

summary(est.m.buildup.bal[[1]])
## 
## ==========================
## Summary of model fit
## ==========================
## 
## Formula:   nw ~ edges + degrange(from = 2) + offset(nodematch("role.class", 
##     diff = TRUE, keep = 1:2))
## <environment: 0x560796482160>
## 
## Iterations:  163 out of 400 
## 
## Monte Carlo MLE Results:
##                        Estimate Std. Error MCMC % p-value    
## edges                  -8.74395    0.03417      0  <1e-04 ***
## deg2+                      -Inf    0.00000      0  <1e-04 ***
## nodematch.role.class.I     -Inf    0.00000      0  <1e-04 ***
## nodematch.role.class.R     -Inf    0.00000      0  <1e-04 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log-likelihood was not estimated for this fit.
## To get deviances, AIC, and/or BIC from fit `object$fit` run 
##   > object$fit<-logLik(object$fit, add=TRUE)
## to add it to the object or rerun this function with eval.loglik=TRUE.
## 
##  Warning: The following terms have infinite coefficient estimates:
##   deg2+ 
## 
##  The following terms are fixed by offset and are not estimated:
##   nodematch.role.class.I nodematch.role.class.R 
## 
## 
## Dissolution Coefficients
## =======================
## Dissolution Model: ~offset(edges)
## Target Statistics: 153
## Crude Coefficient: 5.023881
## Mortality/Exit Rate: 5.302239e-05
## Adjusted Coefficient: 5.040238

Model 2

summary(est.m.buildup.bal[[2]])
## 
## ==========================
## Summary of model fit
## ==========================
## 
## Formula:   nw ~ edges + nodefactor("race..wa", base = 3) + degrange(from = 2) + 
##     offset(nodematch("role.class", diff = TRUE, keep = 1:2))
## <environment: 0x5607ad31cec8>
## 
## Iterations:  102 out of 400 
## 
## Monte Carlo MLE Results:
##                        Estimate Std. Error MCMC % p-value    
## edges                  -8.78864    0.03912      0  <1e-04 ***
## nodefactor.race..wa.B  -0.38471    0.08766      0  <1e-04 ***
## nodefactor.race..wa.H   0.41133    0.06572      0  <1e-04 ***
## deg2+                      -Inf    0.00000      0  <1e-04 ***
## nodematch.role.class.I     -Inf    0.00000      0  <1e-04 ***
## nodematch.role.class.R     -Inf    0.00000      0  <1e-04 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log-likelihood was not estimated for this fit.
## To get deviances, AIC, and/or BIC from fit `object$fit` run 
##   > object$fit<-logLik(object$fit, add=TRUE)
## to add it to the object or rerun this function with eval.loglik=TRUE.
## 
##  Warning: The following terms have infinite coefficient estimates:
##   deg2+ 
## 
##  The following terms are fixed by offset and are not estimated:
##   nodematch.role.class.I nodematch.role.class.R 
## 
## 
## Dissolution Coefficients
## =======================
## Dissolution Model: ~offset(edges)
## Target Statistics: 153
## Crude Coefficient: 5.023881
## Mortality/Exit Rate: 5.302239e-05
## Adjusted Coefficient: 5.040238

Model 3

summary(est.m.buildup.bal[[3]])
## 
## ==========================
## Summary of model fit
## ==========================
## 
## Formula:   nw ~ edges + nodefactor("race..wa", base = 3) + nodematch("race..wa", 
##     diff = TRUE) + degrange(from = 2) + offset(nodematch("role.class", 
##     diff = TRUE, keep = 1:2))
## <environment: 0x5607c1d26e88>
## 
## Iterations:  235 out of 400 
## 
## Monte Carlo MLE Results:
##                        Estimate Std. Error MCMC %  p-value    
## edges                  -10.2727     0.1886      0  < 1e-04 ***
## nodefactor.race..wa.B    0.5925     0.1682      0 0.000426 ***
## nodefactor.race..wa.H    1.4058     0.1830      0  < 1e-04 ***
## nodematch.race..wa.B     1.3391     0.2618      0  < 1e-04 ***
## nodematch.race..wa.H     0.6582     0.2076      0 0.001525 ** 
## nodematch.race..wa.O     1.5656     0.1882      0  < 1e-04 ***
## deg2+                      -Inf     0.0000      0  < 1e-04 ***
## nodematch.role.class.I     -Inf     0.0000      0  < 1e-04 ***
## nodematch.role.class.R     -Inf     0.0000      0  < 1e-04 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log-likelihood was not estimated for this fit.
## To get deviances, AIC, and/or BIC from fit `object$fit` run 
##   > object$fit<-logLik(object$fit, add=TRUE)
## to add it to the object or rerun this function with eval.loglik=TRUE.
## 
##  Warning: The following terms have infinite coefficient estimates:
##   deg2+ 
## 
##  The following terms are fixed by offset and are not estimated:
##   nodematch.role.class.I nodematch.role.class.R 
## 
## 
## Dissolution Coefficients
## =======================
## Dissolution Model: ~offset(edges)
## Target Statistics: 153
## Crude Coefficient: 5.023881
## Mortality/Exit Rate: 5.302239e-05
## Adjusted Coefficient: 5.040238

Model 4

summary(est.m.buildup.bal[[4]])
## 
## ==========================
## Summary of model fit
## ==========================
## 
## Formula:   nw ~ edges + nodefactor("deg.pers") + nodefactor("race..wa", 
##     base = 3) + nodematch("race..wa", diff = TRUE) + degrange(from = 2) + 
##     offset(nodematch("role.class", diff = TRUE, keep = 1:2))
## <environment: 0x5607d6956978>
## 
## Iterations:  152 out of 400 
## 
## Monte Carlo MLE Results:
##                         Estimate Std. Error MCMC % p-value    
## edges                  -10.17034    0.19083      0 < 1e-04 ***
## nodefactor.deg.pers.1   -0.30470    0.06226      0 < 1e-04 ***
## nodefactor.deg.pers.2   -0.09834    0.05912      0 0.09622 .  
## nodefactor.race..wa.B    0.59675    0.16792      0 0.00038 ***
## nodefactor.race..wa.H    1.41065    0.18424      0 < 1e-04 ***
## nodematch.race..wa.B     1.33978    0.25705      0 < 1e-04 ***
## nodematch.race..wa.H     0.66351    0.20730      0 0.00137 ** 
## nodematch.race..wa.O     1.56433    0.18863      0 < 1e-04 ***
## deg2+                       -Inf    0.00000      0 < 1e-04 ***
## nodematch.role.class.I      -Inf    0.00000      0 < 1e-04 ***
## nodematch.role.class.R      -Inf    0.00000      0 < 1e-04 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log-likelihood was not estimated for this fit.
## To get deviances, AIC, and/or BIC from fit `object$fit` run 
##   > object$fit<-logLik(object$fit, add=TRUE)
## to add it to the object or rerun this function with eval.loglik=TRUE.
## 
##  Warning: The following terms have infinite coefficient estimates:
##   deg2+ 
## 
##  The following terms are fixed by offset and are not estimated:
##   nodematch.role.class.I nodematch.role.class.R 
## 
## 
## Dissolution Coefficients
## =======================
## Dissolution Model: ~offset(edges)
## Target Statistics: 153
## Crude Coefficient: 5.023881
## Mortality/Exit Rate: 5.302239e-05
## Adjusted Coefficient: 5.040238

Model 5

summary(est.m.buildup.bal[[5]])
## 
## ==========================
## Summary of model fit
## ==========================
## 
## Formula:   nw ~ edges + nodefactor("deg.pers") + nodefactor("race..wa", 
##     base = 3) + nodefactor("region", base = 2) + nodematch("race..wa", 
##     diff = TRUE) + degrange(from = 2) + offset(nodematch("role.class", 
##     diff = TRUE, keep = 1:2))
## <environment: 0x5607eb6a3ed8>
## 
## Iterations:  106 out of 400 
## 
## Monte Carlo MLE Results:
##                        Estimate Std. Error MCMC %  p-value    
## edges                  -9.86114    0.19286      0  < 1e-04 ***
## nodefactor.deg.pers.1  -0.31120    0.06253      0  < 1e-04 ***
## nodefactor.deg.pers.2  -0.10049    0.05977      0 0.092725 .  
## nodefactor.race..wa.B   0.55488    0.16763      0 0.000932 ***
## nodefactor.race..wa.H   1.42475    0.18313      0  < 1e-04 ***
## nodefactor.region.EW   -0.25065    0.06989      0 0.000335 ***
## nodefactor.region.OW   -0.39516    0.04470      0  < 1e-04 ***
## nodematch.race..wa.B    1.34037    0.25998      0  < 1e-04 ***
## nodematch.race..wa.H    0.66006    0.20699      0 0.001428 ** 
## nodematch.race..wa.O    1.56595    0.18757      0  < 1e-04 ***
## deg2+                      -Inf    0.00000      0  < 1e-04 ***
## nodematch.role.class.I     -Inf    0.00000      0  < 1e-04 ***
## nodematch.role.class.R     -Inf    0.00000      0  < 1e-04 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log-likelihood was not estimated for this fit.
## To get deviances, AIC, and/or BIC from fit `object$fit` run 
##   > object$fit<-logLik(object$fit, add=TRUE)
## to add it to the object or rerun this function with eval.loglik=TRUE.
## 
##  Warning: The following terms have infinite coefficient estimates:
##   deg2+ 
## 
##  The following terms are fixed by offset and are not estimated:
##   nodematch.role.class.I nodematch.role.class.R 
## 
## 
## Dissolution Coefficients
## =======================
## Dissolution Model: ~offset(edges)
## Target Statistics: 153
## Crude Coefficient: 5.023881
## Mortality/Exit Rate: 5.302239e-05
## Adjusted Coefficient: 5.040238

Model 6

summary(est.m.buildup.bal[[6]])
## 
## ==========================
## Summary of model fit
## ==========================
## 
## Formula:   nw ~ edges + nodefactor("deg.pers") + nodefactor("race..wa", 
##     base = 3) + nodefactor("region", base = 2) + nodematch("race..wa", 
##     diff = TRUE) + absdiff("sqrt.age") + degrange(from = 2) + 
##     offset(nodematch("role.class", diff = TRUE, keep = 1:2))
## <environment: 0x560800505398>
## 
## Iterations:  76 out of 400 
## 
## Monte Carlo MLE Results:
##                        Estimate Std. Error MCMC % p-value    
## edges                  -8.73835    0.19472      0 < 1e-04 ***
## nodefactor.deg.pers.1  -0.31083    0.06328      0 < 1e-04 ***
## nodefactor.deg.pers.2  -0.10403    0.05972      0 0.08149 .  
## nodefactor.race..wa.B   0.54964    0.16727      0 0.00102 ** 
## nodefactor.race..wa.H   1.42371    0.18384      0 < 1e-04 ***
## nodefactor.region.EW   -0.24138    0.06966      0 0.00053 ***
## nodefactor.region.OW   -0.39362    0.04476      0 < 1e-04 ***
## nodematch.race..wa.B    1.34183    0.25914      0 < 1e-04 ***
## nodematch.race..wa.H    0.66048    0.20759      0 0.00146 ** 
## nodematch.race..wa.O    1.56524    0.18773      0 < 1e-04 ***
## absdiff.sqrt.age       -1.39760    0.04187      0 < 1e-04 ***
## deg2+                      -Inf    0.00000      0 < 1e-04 ***
## nodematch.role.class.I     -Inf    0.00000      0 < 1e-04 ***
## nodematch.role.class.R     -Inf    0.00000      0 < 1e-04 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log-likelihood was not estimated for this fit.
## To get deviances, AIC, and/or BIC from fit `object$fit` run 
##   > object$fit<-logLik(object$fit, add=TRUE)
## to add it to the object or rerun this function with eval.loglik=TRUE.
## 
##  Warning: The following terms have infinite coefficient estimates:
##   deg2+ 
## 
##  The following terms are fixed by offset and are not estimated:
##   nodematch.role.class.I nodematch.role.class.R 
## 
## 
## Dissolution Coefficients
## =======================
## Dissolution Model: ~offset(edges)
## Target Statistics: 153
## Crude Coefficient: 5.023881
## Mortality/Exit Rate: 5.302239e-05
## Adjusted Coefficient: 5.040238

Model 7

summary(est.m.buildup.bal[[7]])
## 
## ==========================
## Summary of model fit
## ==========================
## 
## Formula:   nw ~ edges + nodefactor("deg.pers") + nodefactor("race..wa", 
##     base = 3) + nodefactor("region", base = 2) + nodematch("race..wa", 
##     diff = TRUE) + absdiff("sqrt.age") + degrange(from = 2) + 
##     offset(nodematch("role.class", diff = TRUE, keep = 1:2)) + 
##     offset(nodemix("region", base = c(1, 3, 6)))
## <environment: 0x560815414bc8>
## 
## Iterations:  59 out of 400 
## 
## Monte Carlo MLE Results:
##                         Estimate Std. Error MCMC %  p-value    
## edges                  -8.315317   0.194012      0  < 1e-04 ***
## nodefactor.deg.pers.1  -0.310127   0.062656      0  < 1e-04 ***
## nodefactor.deg.pers.2  -0.104764   0.059749      0 0.079530 .  
## nodefactor.race..wa.B   0.606938   0.168067      0 0.000305 ***
## nodefactor.race..wa.H   1.494859   0.183869      0  < 1e-04 ***
## nodefactor.region.EW    0.666595   0.060014      0  < 1e-04 ***
## nodefactor.region.OW   -0.009103   0.037916      0 0.810263    
## nodematch.race..wa.B    1.216911   0.260142      0  < 1e-04 ***
## nodematch.race..wa.H    0.509008   0.208100      0 0.014446 *  
## nodematch.race..wa.O    1.627132   0.188217      0  < 1e-04 ***
## absdiff.sqrt.age       -1.397501   0.041976      0  < 1e-04 ***
## deg2+                       -Inf   0.000000      0  < 1e-04 ***
## nodematch.role.class.I      -Inf   0.000000      0  < 1e-04 ***
## nodematch.role.class.R      -Inf   0.000000      0  < 1e-04 ***
## mix.region.EW.KC            -Inf   0.000000      0  < 1e-04 ***
## mix.region.EW.OW            -Inf   0.000000      0  < 1e-04 ***
## mix.region.KC.OW            -Inf   0.000000      0  < 1e-04 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log-likelihood was not estimated for this fit.
## To get deviances, AIC, and/or BIC from fit `object$fit` run 
##   > object$fit<-logLik(object$fit, add=TRUE)
## to add it to the object or rerun this function with eval.loglik=TRUE.
## 
##  Warning: The following terms have infinite coefficient estimates:
##   deg2+ 
## 
##  The following terms are fixed by offset and are not estimated:
##   nodematch.role.class.I nodematch.role.class.R mix.region.EW.KC mix.region.EW.OW mix.region.KC.OW 
## 
## 
## Dissolution Coefficients
## =======================
## Dissolution Model: ~offset(edges)
## Target Statistics: 153
## Crude Coefficient: 5.023881
## Mortality/Exit Rate: 5.302239e-05
## Adjusted Coefficient: 5.040238

Model 8

summary(est.m.buildup.bal[[8]])
## 
## ==========================
## Summary of model fit
## ==========================
## 
## Formula:   nw ~ edges + nodefactor("deg.pers") + nodefactor("race..wa", 
##     base = 3) + nodefactor("region", base = 2) + nodematch("race..wa", 
##     diff = TRUE) + absdiff("sqrt.age") + nodematch("region", 
##     diff = FALSE) + degrange(from = 2) + offset(nodematch("role.class", 
##     diff = TRUE, keep = 1:2))
## <environment: 0x56082a4104c0>
## 
## Iterations:  87 out of 400 
## 
## Monte Carlo MLE Results:
##                         Estimate Std. Error MCMC %  p-value    
## edges                  -11.03239    0.20677      0  < 1e-04 ***
## nodefactor.deg.pers.1   -0.31090    0.06278      0  < 1e-04 ***
## nodefactor.deg.pers.2   -0.10419    0.05970      0 0.080933 .  
## nodefactor.race..wa.B    0.59298    0.16847      0 0.000432 ***
## nodefactor.race..wa.H    1.47743    0.18489      0  < 1e-04 ***
## nodefactor.region.EW     0.56263    0.06247      0  < 1e-04 ***
## nodefactor.region.OW    -0.04204    0.03881      0 0.278705    
## nodematch.race..wa.B     1.24035    0.26056      0  < 1e-04 ***
## nodematch.race..wa.H     0.54183    0.20893      0 0.009503 ** 
## nodematch.race..wa.O     1.61170    0.18837      0  < 1e-04 ***
## absdiff.sqrt.age        -1.39797    0.04217      0  < 1e-04 ***
## nodematch.region         2.66475    0.07157      0  < 1e-04 ***
## deg2+                       -Inf    0.00000      0  < 1e-04 ***
## nodematch.role.class.I      -Inf    0.00000      0  < 1e-04 ***
## nodematch.role.class.R      -Inf    0.00000      0  < 1e-04 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log-likelihood was not estimated for this fit.
## To get deviances, AIC, and/or BIC from fit `object$fit` run 
##   > object$fit<-logLik(object$fit, add=TRUE)
## to add it to the object or rerun this function with eval.loglik=TRUE.
## 
##  Warning: The following terms have infinite coefficient estimates:
##   deg2+ 
## 
##  The following terms are fixed by offset and are not estimated:
##   nodematch.role.class.I nodematch.role.class.R 
## 
## 
## Dissolution Coefficients
## =======================
## Dissolution Model: ~offset(edges)
## Target Statistics: 153
## Crude Coefficient: 5.023881
## Mortality/Exit Rate: 5.302239e-05
## Adjusted Coefficient: 5.040238

Network diagnostics

Model 1

(dx_main1 <- netdx(est.m.buildup.bal[[1]], nsims = 10, nsteps = 1000, ncores = 4, nwstats.formula = est.m.buildup.bal[[8]]$formation))
## 
## Network Diagnostics
## -----------------------
## - Simulating 10 networks
## - Calculating formation statistics
## - Calculating duration statistics
## - Calculating dissolution statistics
## 
## EpiModel Network Diagnostics
## =======================
## Diagnostic Method: Dynamic
## Simulations: 10
## Time Steps per Sim: 1000
## 
## Formation Diagnostics
## ----------------------- 
##                        Target Sim Mean Pct Diff Sim SD
## edges                  2240.5 2234.827   -0.003 31.719
## nodefactor.deg.pers.1      NA  572.696       NA 17.693
## nodefactor.deg.pers.2      NA  616.343       NA 18.859
## nodefactor.race..wa.B      NA  272.185       NA 12.206
## nodefactor.race..wa.H      NA  483.940       NA 16.751
## nodefactor.region.EW       NA  453.268       NA 16.519
## nodefactor.region.OW       NA 1462.782       NA 30.578
## nodematch.race..wa.B       NA    8.232       NA  2.685
## nodematch.race..wa.H       NA   26.348       NA  4.634
## nodematch.race..wa.O       NA 1543.153       NA 26.422
## absdiff.sqrt.age           NA 2547.696       NA 51.862
## nodematch.region           NA  991.177       NA 26.059
## deg2+                      NA    0.000       NA  0.000
## nodematch.role.class.I     NA    0.000       NA  0.000
## nodematch.role.class.R     NA    0.000       NA  0.000
## 
## Dissolution Diagnostics
## ----------------------- 
##                 Target Sim Mean Pct Diff  Sim SD
## Edge Duration  153.000  125.703   -0.178 120.246
## Pct Edges Diss   0.007    0.007    0.003   0.002
plot(dx_main1, type="formation")

plot(dx_main1, type="duration")

plot(dx_main1, type="dissolution")

Model 2

(dx_main2 <- netdx(est.m.buildup.bal[[2]], nsims = 10, nsteps = 1000, ncores = 4, nwstats.formula = est.m.buildup.bal[[8]]$formation))
## 
## Network Diagnostics
## -----------------------
## - Simulating 10 networks
## - Calculating formation statistics
## - Calculating duration statistics
## - Calculating dissolution statistics
## 
## EpiModel Network Diagnostics
## =======================
## Diagnostic Method: Dynamic
## Simulations: 10
## Time Steps per Sim: 1000
## 
## Formation Diagnostics
## ----------------------- 
##                          Target Sim Mean Pct Diff Sim SD
## edges                  2240.500 2242.452    0.001 28.395
## nodefactor.deg.pers.1        NA  574.868       NA 17.992
## nodefactor.deg.pers.2        NA  619.263       NA 18.074
## nodefactor.race..wa.B   213.834  216.734    0.014 12.680
## nodefactor.race..wa.H   587.844  589.952    0.004 15.627
## nodefactor.region.EW         NA  469.565       NA 17.691
## nodefactor.region.OW         NA 1469.677       NA 30.724
## nodematch.race..wa.B         NA    5.244       NA  1.831
## nodematch.race..wa.H         NA   36.657       NA  6.049
## nodematch.race..wa.O         NA 1505.443       NA 25.217
## absdiff.sqrt.age             NA 2548.298       NA 53.578
## nodematch.region             NA  986.087       NA 26.643
## deg2+                        NA    0.000       NA  0.000
## nodematch.role.class.I       NA    0.000       NA  0.000
## nodematch.role.class.R       NA    0.000       NA  0.000
## 
## Dissolution Diagnostics
## ----------------------- 
##                 Target Sim Mean Pct Diff  Sim SD
## Edge Duration  153.000  126.149   -0.175 120.750
## Pct Edges Diss   0.007    0.007   -0.003   0.002
plot(dx_main2, type="formation")

plot(dx_main2, type="duration")

plot(dx_main2, type="dissolution")

Model 3

(dx_main3 <- netdx(est.m.buildup.bal[[3]], nsims = 10, nsteps = 1000, ncores = 4, nwstats.formula = est.m.buildup.bal[[8]]$formation))
## 
## Network Diagnostics
## -----------------------
## - Simulating 10 networks
## - Calculating formation statistics
## - Calculating duration statistics
## - Calculating dissolution statistics
## 
## EpiModel Network Diagnostics
## =======================
## Diagnostic Method: Dynamic
## Simulations: 10
## Time Steps per Sim: 1000
## 
## Formation Diagnostics
## ----------------------- 
##                          Target Sim Mean Pct Diff Sim SD
## edges                  2240.500 2237.948   -0.001 30.837
## nodefactor.deg.pers.1        NA  569.327       NA 16.739
## nodefactor.deg.pers.2        NA  616.123       NA 19.700
## nodefactor.race..wa.B   213.834  212.875   -0.004 12.682
## nodefactor.race..wa.H   587.844  580.438   -0.013 17.648
## nodefactor.region.EW         NA  466.720       NA 15.794
## nodefactor.region.OW         NA 1468.944       NA 27.578
## nodematch.race..wa.B     31.177   31.688    0.016  5.116
## nodematch.race..wa.H    123.300  118.437   -0.039  9.352
## nodematch.race..wa.O   1638.946 1640.231    0.001 26.542
## absdiff.sqrt.age             NA 2538.402       NA 49.974
## nodematch.region             NA  986.999       NA 28.130
## deg2+                        NA    0.000       NA  0.000
## nodematch.role.class.I       NA    0.000       NA  0.000
## nodematch.role.class.R       NA    0.000       NA  0.000
## 
## Dissolution Diagnostics
## ----------------------- 
##                 Target Sim Mean Pct Diff  Sim SD
## Edge Duration  153.000  126.239   -0.175 121.093
## Pct Edges Diss   0.007    0.007   -0.003   0.002
plot(dx_main3, type="formation")

plot(dx_main3, type="duration")

plot(dx_main3, type="dissolution")

Model 4

(dx_main4 <- netdx(est.m.buildup.bal[[4]], nsims = 10, nsteps = 1000, ncores = 4, nwstats.formula = est.m.buildup.bal[[8]]$formation))
## 
## Network Diagnostics
## -----------------------
## - Simulating 10 networks
## - Calculating formation statistics
## - Calculating duration statistics
## - Calculating dissolution statistics
## 
## EpiModel Network Diagnostics
## =======================
## Diagnostic Method: Dynamic
## Simulations: 10
## Time Steps per Sim: 1000
## 
## Formation Diagnostics
## ----------------------- 
##                          Target Sim Mean Pct Diff Sim SD
## edges                  2240.500 2234.744   -0.003 28.619
## nodefactor.deg.pers.1   493.000  494.521    0.003 18.906
## nodefactor.deg.pers.2   603.000  602.727    0.000 19.572
## nodefactor.race..wa.B   213.834  208.750   -0.024 12.620
## nodefactor.race..wa.H   587.844  580.989   -0.012 17.654
## nodefactor.region.EW         NA  466.656       NA 16.416
## nodefactor.region.OW         NA 1467.211       NA 30.367
## nodematch.race..wa.B     31.177   30.292   -0.028  5.077
## nodematch.race..wa.H    123.300  119.041   -0.035  8.643
## nodematch.race..wa.O   1638.946 1639.353    0.000 25.838
## absdiff.sqrt.age             NA 2536.368       NA 45.469
## nodematch.region             NA  983.803       NA 25.381
## deg2+                        NA    0.000       NA  0.000
## nodematch.role.class.I       NA    0.000       NA  0.000
## nodematch.role.class.R       NA    0.000       NA  0.000
## 
## Dissolution Diagnostics
## ----------------------- 
##                 Target Sim Mean Pct Diff  Sim SD
## Edge Duration  153.000  125.719   -0.178 120.412
## Pct Edges Diss   0.007    0.007    0.000   0.002
plot(dx_main4, type="formation")

plot(dx_main4, type="duration")

plot(dx_main4, type="dissolution")

Model 5

(dx_main5 <- netdx(est.m.buildup.bal[[5]], nsims = 10, nsteps = 1000, ncores = 4, nwstats.formula = est.m.buildup.bal[[8]]$formation))
## 
## Network Diagnostics
## -----------------------
## - Simulating 10 networks
## - Calculating formation statistics
## - Calculating duration statistics
## - Calculating dissolution statistics
## 
## EpiModel Network Diagnostics
## =======================
## Diagnostic Method: Dynamic
## Simulations: 10
## Time Steps per Sim: 1000
## 
## Formation Diagnostics
## ----------------------- 
##                          Target Sim Mean Pct Diff Sim SD
## edges                  2240.500 2235.535   -0.002 28.357
## nodefactor.deg.pers.1   493.000  493.286    0.001 16.312
## nodefactor.deg.pers.2   603.000  600.457   -0.004 16.428
## nodefactor.race..wa.B   213.834  212.368   -0.007 12.291
## nodefactor.race..wa.H   587.844  579.359   -0.014 16.352
## nodefactor.region.EW    445.561  448.030    0.006 14.550
## nodefactor.region.OW   1278.131 1277.676    0.000 28.187
## nodematch.race..wa.B     31.177   31.341    0.005  4.909
## nodematch.race..wa.H    123.300  117.559   -0.047  9.629
## nodematch.race..wa.O   1638.946 1637.982   -0.001 26.233
## absdiff.sqrt.age             NA 2546.125       NA 46.089
## nodematch.region             NA 1046.027       NA 21.525
## deg2+                        NA    0.000       NA  0.000
## nodematch.role.class.I       NA    0.000       NA  0.000
## nodematch.role.class.R       NA    0.000       NA  0.000
## 
## Dissolution Diagnostics
## ----------------------- 
##                 Target Sim Mean Pct Diff  Sim SD
## Edge Duration  153.000  125.548   -0.179 120.292
## Pct Edges Diss   0.007    0.007    0.002   0.002
plot(dx_main5, type="formation")

plot(dx_main5, type="duration")

plot(dx_main5, type="dissolution")

Model 6

(dx_main6 <- netdx(est.m.buildup.bal[[6]], nsims = 10, nsteps = 1000, ncores = 4, nwstats.formula = est.m.buildup.bal[[8]]$formation))
## 
## Network Diagnostics
## -----------------------
## - Simulating 10 networks
## - Calculating formation statistics
## - Calculating duration statistics
## - Calculating dissolution statistics
## 
## EpiModel Network Diagnostics
## =======================
## Diagnostic Method: Dynamic
## Simulations: 10
## Time Steps per Sim: 1000
## 
## Formation Diagnostics
## ----------------------- 
##                          Target Sim Mean Pct Diff Sim SD
## edges                  2240.500 2238.102   -0.001 28.109
## nodefactor.deg.pers.1   493.000  490.743   -0.005 16.184
## nodefactor.deg.pers.2   603.000  601.905   -0.002 18.063
## nodefactor.race..wa.B   213.834  214.122    0.001 12.458
## nodefactor.race..wa.H   587.844  573.825   -0.024 16.327
## nodefactor.region.EW    445.561  438.963   -0.015 14.526
## nodefactor.region.OW   1278.131 1283.462    0.004 28.911
## nodematch.race..wa.B     31.177   31.260    0.003  4.311
## nodematch.race..wa.H    123.300  112.323   -0.089  7.582
## nodematch.race..wa.O   1638.946 1639.690    0.000 25.723
## absdiff.sqrt.age       1206.285 1212.373    0.005 30.176
## nodematch.region             NA 1059.746       NA 22.519
## deg2+                        NA    0.000       NA  0.000
## nodematch.role.class.I       NA    0.000       NA  0.000
## nodematch.role.class.R       NA    0.000       NA  0.000
## 
## Dissolution Diagnostics
## ----------------------- 
##                 Target Sim Mean Pct Diff  Sim SD
## Edge Duration  153.000  126.405   -0.174 120.702
## Pct Edges Diss   0.007    0.007   -0.004   0.002
plot(dx_main6, type="formation")

plot(dx_main6, type="duration")

plot(dx_main6, type="dissolution")

Model 7

(dx_main7 <- netdx(est.m.buildup.bal[[7]], nsims = 10, nsteps = 1000, ncores = 4))
## 
## Network Diagnostics
## -----------------------
## - Simulating 10 networks
## - Calculating formation statistics
## - Calculating duration statistics
## - Calculating dissolution statistics
## 
## EpiModel Network Diagnostics
## =======================
## Diagnostic Method: Dynamic
## Simulations: 10
## Time Steps per Sim: 1000
## 
## Formation Diagnostics
## ----------------------- 
##                          Target Sim Mean Pct Diff Sim SD
## edges                  2240.500 2205.330   -0.016 27.924
## nodefactor.deg.pers.1   493.000  489.719   -0.007 16.360
## nodefactor.deg.pers.2   603.000  593.763   -0.015 17.924
## nodefactor.race..wa.B   213.834  212.831   -0.005 12.175
## nodefactor.race..wa.H   587.844  543.241   -0.076 17.768
## nodefactor.region.EW    445.561  408.685   -0.083 18.079
## nodefactor.region.OW   1278.131 1266.092   -0.009 34.580
## nodematch.race..wa.B     31.177   30.025   -0.037  4.719
## nodematch.race..wa.H    123.300   92.966   -0.246  8.569
## nodematch.race..wa.O   1638.946 1617.693   -0.013 26.245
## absdiff.sqrt.age       1206.285 1213.921    0.006 25.013
## deg2+                        NA    0.000       NA  0.000
## nodematch.role.class.I       NA    0.000       NA  0.000
## nodematch.role.class.R       NA    0.000       NA  0.000
## mix.region.EW.KC             NA    0.000       NA  0.000
## mix.region.EW.OW             NA    0.000       NA  0.000
## mix.region.KC.OW             NA    0.000       NA  0.000
## 
## Dissolution Diagnostics
## ----------------------- 
##                 Target Sim Mean Pct Diff  Sim SD
## Edge Duration  153.000  125.720   -0.178 120.462
## Pct Edges Diss   0.007    0.007   -0.003   0.002
plot(dx_main7, type="formation")

plot(dx_main7, type="duration")

plot(dx_main7, type="dissolution")

Model 8

(dx_main8 <- netdx(est.m.buildup.bal[[7]], nsims = 10, nsteps = 1000, ncores = 4))
## 
## Network Diagnostics
## -----------------------
## - Simulating 10 networks
## - Calculating formation statistics
## - Calculating duration statistics
## - Calculating dissolution statistics
## 
## EpiModel Network Diagnostics
## =======================
## Diagnostic Method: Dynamic
## Simulations: 10
## Time Steps per Sim: 1000
## 
## Formation Diagnostics
## ----------------------- 
##                          Target Sim Mean Pct Diff Sim SD
## edges                  2240.500 2202.823   -0.017 30.321
## nodefactor.deg.pers.1   493.000  491.113   -0.004 18.044
## nodefactor.deg.pers.2   603.000  593.405   -0.016 17.964
## nodefactor.race..wa.B   213.834  210.918   -0.014 11.593
## nodefactor.race..wa.H   587.844  547.518   -0.069 18.066
## nodefactor.region.EW    445.561  404.939   -0.091 18.431
## nodefactor.region.OW   1278.131 1262.948   -0.012 30.986
## nodematch.race..wa.B     31.177   29.479   -0.054  4.739
## nodematch.race..wa.H    123.300   94.298   -0.235  7.839
## nodematch.race..wa.O   1638.946 1615.399   -0.014 26.837
## absdiff.sqrt.age       1206.285 1211.555    0.004 28.374
## deg2+                        NA    0.000       NA  0.000
## nodematch.role.class.I       NA    0.000       NA  0.000
## nodematch.role.class.R       NA    0.000       NA  0.000
## mix.region.EW.KC             NA    0.000       NA  0.000
## mix.region.EW.OW             NA    0.000       NA  0.000
## mix.region.KC.OW             NA    0.000       NA  0.000
## 
## Dissolution Diagnostics
## ----------------------- 
##                 Target Sim Mean Pct Diff  Sim SD
## Edge Duration  153.000  125.410   -0.180 120.179
## Pct Edges Diss   0.007    0.007   -0.001   0.002
plot(dx_main8, type="formation")

plot(dx_main8, type="duration")

plot(dx_main8, type="dissolution")